Revista 

CIENCIA 



Vol. 13, 1, 1-18 (2010) 



Revista 

CIENCIA 

EDITOR 

Roberto Aguiar Falcom 

Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador 
e-mail: raguiar® espe.edu.ee 

COMITE EDITORIAL 



A.H. Barbat 

E.T.S. Ingenieros de Caminos y Canales y Puertos 
Universidad Politecnica de Cantuna 
Campus Norte UPC, 08034. 


A. Tierra 

Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador. 


J. L. Almazan 

Escuela de Ingenieria 

Pontificia Universidad Catolica de Santiago 
Santiago de Chile 


M- Lafuente 

Instituto de Materiales y Modelos Estructurales 
Universidad Central de Venezuela 
Caracas, Venezuela. 


L. Cordova 

Facultad de Ciencias y Humanidades 
Universidad del Valle de Guatemala 
Guatemala. 


W. Pozo 

Centro de Investigaciones del lASA 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador. 


P. Dechent 

Escuela de Ingenieria 
Universidad de Concepcion 
Concepcion, Chile 


M. Romo 

Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador. 


S. Fernandez 

Academia de Estructuras 
Instituto Politecnico Nacional 
Mexico D.F. 


J. M Ruiz 

Facultad de Construcciones 
Departamento de Ingenieria Civil 
Universidad de Oriente 
Santiago de Cuba, Cuba. 


P. Caiza 

Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador 


S. Ruiz 

Instituto de Ingenieria 

Universidad Nacional Autonoma de Mexico 
Mexico D.F. 


C. Genatios 

Instituto de Materiales y Modelos Estructurales 
Universidad Central de Venezuela 
Caracas, Venezuela. 


F. Williams L. 

Instituto de Ingenieria 
Universidad Veracruzana 
Veracruz, Mexico. 


L. Cumbal 

Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador. 

FORMATO DE TEXTO 

Javier Villamarm 

Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador 


A. Oleas 

Centro de Investigaciones del lASA 
Escuela Politecnica del Ejercito 
Valle de los Chillos, Ecuador. 



ISSN 1390-1117 

©2010 ESPE, Quito, Ecuador 



2 



Revista CIENCIA Volumen 13 N.- 1 de 2010 



Revista 

CIENCIA 

Sumario 

Volumen 13, mimero 1, 2010 

Leyes de Atenuacion para Sismos Corticales y de Subduccion para el Ecuador 001 

Roberto Aguiar Falconi, Edwin Garcia y Javier Villamarin 



Analisis Sismico en sentido longitudinal del Puente Sur que tiene aisladores EPS con 
CEINCI-LAB. Comparacion entre dos modelos 019 

Roberto Aguiar Falconi y Darwin Carrera 



Pesticidas Naturales y Sinteticos 043 

Vladimir Aguirre Yela y Vicente Delgado 



Evaluacion de dos metodos para la reduccion de nitrogeno, fosforo y DQO de aguas 
residuales, mediante un cultivo axenico de cianobacterias y un consorcio microbiano, 
inmovilizados y en suspension 055 

Mauricio Moreno, Blanca Naranjo y Alma Koch 



Analisis sismico en sentido longitudinal del Puente Carrizal con dos modelos. 

Empleando CEINCI-LAB 063 

Roberto Aguiar Falconi y Edwin Garcia 



Utilizacion de bongos micorricico-arbusculares en plantulas de tomate de arbol 
(solanum betaceum) cultivadas in vitrei- Efectos durante la fase de aclimatacion 087 

Espin E., Medina M.E., Jadan M, & Proano K. 



Evaluacion de la actividad del formononetin, acido 3-indol buritico y tiamina 
para promover la germinacion in vitro de esporas de bongos 

micorricico-arbusculares (HMA) 095 

Baroja 1., Medina M.E. & Proano K. 



Revista Semestral de la Escuela Politecnica del Ejercito Quito, Ecuador. 
ISSN 1390-1117 
© 2010 ESPE, Quito — Ecuador 
ESPE, Quito — Ecuador 




Revista CIENCIA Volumen 13 N.- 1 de 2010 



3 



LEYES DE ATENUACION PARA SISMOS CORTICALES Y 
DE SUBDUCCION PARA EL ECUADOR 



Roberto Aguiar Falconi, Edwin Garcia y Javier Viiiamarm 



Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
raauiar@espe.edu. ec 



RESUMEN 

En base a la informacion macrosismica de 14 sismos asociados a fallas corticales y a 
7 sismos de subduccion, la mayor parte de ellos registrados en el Ecuador y unos pocos en el 
Sur de Colombia pero que afectaron al Ecuador, se determine una Ley de Atenuacion de 
Intensidades en la escala de Mercalli Modificada. Se tienen 69 dates para los sismos corticales 
y 38 dates para los sismos de subduccion. 

Por otro lade, en base a informacion macrosismica y de registros sismicos, se obtuvo 
una relacion entre la Intensidad maxima en la Escala de Mercalli y la aceleracion maxima del 
suelo, con los dates de 7 sismos registrados en America del Sur. En total se trabajo con 99 
puntos. 



Posteriormente, se presentan dos propuestas de Leyes de Atenuacion de las 
aceleraciones del suelo, para ser utilizados en estudios de Peligrosidad Sismica en el Ecuador. 
La primera encontrada a partir de la relacion entre la aceleracion maxima del suelo y de la 
intensidad maxima, obtenida en este estudio y la segunda empleando la relacion de Saragoni 
et al. (1982). De tal manera que se obtuvieron cuatro Leyes de Atenuacion, dos para la primera 
propuesta y dos para la segunda propuesta. 

Para cada propuesta se tienen dos leyes, una para sismos corticales y la otra para 
sismos de subduccion. Estas propuestas fueron comparadas con otros estudios realizados por 
Aguiar y Castro en 2009; y con leyes obtenidas a partir de registros sismicos en Chile, por Ruiz 
y Saragoni en 2005. 



ABSTRACT 

Based on macro seismic information for 14 earthquakes associated with crustal faults 
and 7 subduction earthquakes, most of them registered in Ecuador and a few in southern 
Colombia but that affected Ecuador, a Law of Attenuation using the Modified Mercalli Intensity 
Scale was determined. There are 69 data bases for cortical earthquakes and 38 for subduction 
earthquakes. 

On the other hand, based on macroseismic information and seismic records -7 
earthquakes recorded in South America-, we obtained a relationship between the maximum 
Intensity on Mercalli scale and the maximum ground acceleration. In total we worked with 99 
points. 



Subsequently, two proposals of Laws of Attenuation for ground accelerations are 
presented for use in studies of Seismic Hazard in Ecuador. The first one was found from the 
relationship between maximum ground acceleration and the maximum intensity obtained in this 
study and the second using Saragoni et al. (1982) relationship. In such a way four Laws of 
Attenuation were obtained, two for the first proposal and two for the second proposal. 
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For each proposal there are two laws, one for crustal earthquakes and the other for 
subductlon earthquakes. These proposals were compared with other studies conducted by 
Aguiar and Castro In 2009 and with laws obtained from seismic records In Chile by Ruiz and 
SaragonI In 2005. 



1. INTRODUCCION 



Venezuela, Colombia, Ecuador, Peru, Bolivia, Chile y Argentina se han visto afectados, 
muy frecuentemente, por grandes sismos que han dejado grandes perdidas humanas y 
economicas. En la tabla 1 se IndIcan los sismos mas destructivos que se han registrado en 
America del Sur en los ultimos 15 ahos; se tiene un total de 14 sismos lo que en promedio 
significa un sismo catastrofico por aho. 



Tabla 1 Sismos destructivos registrados en los ultimos 15 ahos en America del Sur. 



Sismo 


Fecha 


Hora 


Latitud 


Longitud 


Magnitud 


Profundidad 


intensidad 


Ecuador 


28/03/1996 


18:03 


1.04 S 


78.72 W 


Mw= 


5,9 


15,2 km 


VII 


Peru 


12/11/1996 


17:00 


15.30 S 


76.44 W 


Md= 


6,4 


14,0 km 


VII 


Venezuela 


09/07/1997 


15:24 


10.43 S 


63.49 W 


Ms= 


6,8 


10,0 km 


VIII 


Bolivia 


22/05/1998 


1:48 


17.73 S 


65.43 W 


Mw= 


6,8 


24,0 km 


VIII 


Ecuador 


04/08/1 998 


13:59 


0.55 S 


80.53 W 


Ms= 


7,1 


39,0 km 


VIII 


Colombia 


25/01/1999 


13:19 


4.41 N 


75.72 W 


Ml= 


6,2 


<15 km 


VIII 


Peru 


23/06/2001 


15:33 


16.15 S 


73.40 W 


Mw= 


7,4 


33,0 km 


VII-VIII 


Colombia 


15/11/2004 


9:06 


4.81 N 


77.79 W 


Mw= 


7,2 


< 33 km 


VII * 


Chile 


13/06/2005 


22:44 


19.89 N 


69.12 W 


Mw= 


7,9 


117 km 


VII 


Peru 


25/09/2005 


20:56 


5.68 N 


76.4 W 


Mw= 


7,5 


115,0 km 


VI 


Peru 


1 5/08/2007 


18:40 


13.49 S 


76.85 W 


Mw= 


7,9 


26,0 km 


VII-VIII 


Chile 


14/11/2007 


12:40 


22.18 S 


69.84 W 


Mw= 


7,7 


60,0 km 


VIII 


Colombia 


24/05/2008 


19:20 


4.44 N 


73.67 W 


Mw= 


5,9 


4,0 km 


IV* 


Chile 


27/02/2010 


3:34 


35.85 N 


72.72 W 


Mw= 


8,8 


47,4 km 


IX 



* Escala EMS-98 



De los 14 sismos indicados en la tabla 1, 4 son registrados en Peru, 3 en Chile, 3 en 
Colombia, 2 en Ecuador, 1 en Venezuela y 1 en Bolivia. Considerando la longitud de la costa 
del Ecuador con relacion a la longitud de las costas de Colombia, Peru y Chile, el numero de 
sismos registrados es muy alto. Todo esto refleja la gran peligrosidad sfsmica de todos los 
pafses mencionados. 

En la tabla 2 se presenta el numero de vfctimas de los sismos anotados en la tabla 1 . 
Es impresionante la cantidad de muertos que dejaron tres sismos: el que afecto al eje cafetero 
de Colombia el 25 de enero de 1999; el sismo de Pisco de 2007 en Peru y el sismo de 2010 
que causo tanto daho en la region central de Chile. El numero de muertos de estos tres 
terremotos es de 2624 personas. Cantidad sumamente alta por que los sismos no matan gente, 
son las estructuras mal disehadas las que matan a las personas, son la mala ubicacion en que 
se encuentran esas viviendas que han side construidas a lo mejor sin la aprobacion de las 
autoridades municipales, no de otra manera se explica como se permite construir en la misma 
playa a sabiendas de que son lugares muy vulnerables ante un tsunami. 

En el sismo de Colombia de 1999, los Departamentos de Quindfo y Manizales sufrieron 
mucho daho, una muy buena parte de las construcciones que se vinieron a bajo fueron de 
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adobe y otras no tuvieron diseno sismo resistente las cuales fueron construidos antes de 1984, 
ano en el que se promulgo la primera normativa de diseno sismo resistente en Colombia. 
Aguiar (1999). 

La Ciudad de Pisco, luego del sismo de 2007, parecia una cludad que habia sufrido 
gran bombardeo en una guerra; la mayor parte de sus construcclones eran de adobe y no 
habian sido reforzadas. Cludades como Guadalupe, con un 90% de construcclones de adobe 
fueron danadas totalmente, por eso el gran numero de muertos y herldos. Pero tamblen 
edificlos de hormigon armado 



Tabla 2 Perdidas humanas y economicas en los ultimos 15 afios. 



Sismo 


Muertos 


Heridos 


Perdida Economica 


Ecuador 


1996 


27 


156 


5 millones de US 


Peru 


1996 


17 


1591 


43 millones de US 


Venezuela 


1997 


74 


522 


100 millones de US 


Bolivia 


1998 


105 


315 


15 millones de US 


Ecuador 


1998 


3 


52 


30 millones de US 


Colombia 


1999 


1230 


5300 


1800 millones de US 


Peru 


2001 


102 


1368 




Colombia 


2004 








Chile 


2005 


11 


182 


40 millones de US 


Peru 


2005 


10 


164 




Peru 


2007 


595 


1368 


230 millones de US 


Chile 


2007 


2 


140 




Colombia 


2008 


11 


50 




Chile 


2010 


799 




15000 millones de US 




Figura 1 Estructuras danadas en el sismo del Peru de 2007. 




6 



Revista CIENCIA Volumen 13 N.- 1 de 2010 



En la parte superior izquierda de la figura 1 , se aprecia un edificio con p/so blando que 
colapso en la ciudad de Pisco, a la derecha dentro de un circulo se ve el gran dano en la 
columna de un centre de salud, esta falla se denomina columna corta. En la parte inferior 
izquierda se ve el colapso de una casa de adobe de la ciudad de Guadalupe. Finalmente en la 
figura inferior derecha se aprecia el daho debido al cheque de edificios. Todo esto sucedio en 
el sismo de Peru de 2007. 

El ultimo sismo fuerte fue el registrado el 27 de febrero de 2010 en Chile, nuevamente 
las construcciones de adobe fueron destruidas durante el sismo, lo propio sucedio con 
edificaciones sin mantenimiento pero tambien edificios de hormigon armado colapsaron en la 
ciudad de Concepcion. Un buen numero de personas murio con el Tsunami que se dio despues 
del sismo. 

Todo esto nos hace ver la necesidad de investigar en varies campos: el sismologico, 
para poder realizar estudios de peligrosidad sfsmica con la mejor informacion sfsmica 
disponible; en dispositivos de control pasivo como aisladores de base para minimizar la energfa 
sfsmica que ingresa a una estructura y causa daho; en disipadores de energfa pensando en 
reforzar estructuras, entre otros temas. 

En el Ecuador no se dispone de registros con aceleraciones mayores a 0.2 g., siendo 
“g” la aceleracion de la gravedad. Del ultimo sismo fuerte de Bahfa de Caraquez de 1998, no se 
tiene un acelerograma con valor mayor al ya anotado. En los estudios de peligrosidad sfsmica 
se necesita una ley de atenuacion pero al no tener acelerogramas de sismos fuertes, no queda 
otra opcion que el obtener una ley de atenuacion a partir de informacion macrosfsmica y es lo 
que se realize en el presente trabajo. 



Tabla 3 Valores de D para cada Intensidad de isosistas. Aguiar y Castro (2009). 



# 


SISMOS 


FECHA 


Mw. 


D (Km). 


X 


IX 


VIII 


VII 


VI 


V 


IV 


■■ 


II 


1 


RIOBAMBA 


04-02-1797 


7.70 


25.25 


49.14 


88.69 


136.20 


182.40 


230.90 


276.76 


392.23 




2 


OTAVALO 


16-08-1868 


7.17 


44.79 


81.23 


105.31 


233.45 


271.51 










3 


ESMERALDAS 


31/01/1 906 


8.60 




29.48 


84.76 


153.57 


215.94 


291.26 








4 


PASTO CALLE 


1 5/09/1 944 


5.92 








34.03 


38.47 


45.15 








5 


PELILEO 


05/08/1 948 


6.93 


11.46 


14.80 


20.95 


28.94 


37.50 


53.90 


77.24 






6 


AMBATO 


05/08/1949 


6.37 




77.79 


93.89 


143.76 




334.23 








7 


NAPO 


11/05/1 955 


6.98 






40.06 


66.71 


180.72 










8 


BAHIA 


16/01/1956 


7.30 






32.60 


53.00 


107.26 


214.95 


326.78 






9 


ESMERALDAS 


10/01/1958 


7.80 








101.97 


146.89 


219.59 


284.28 






10 


PASA 


30/07/1 960 


5.97 








22.14 


24.28 


30.99 


55.86 


88.02 




11 


ALAUSf 


08/04/1 961 


6.35 








26.63 


43.15 


57.47 


83.58 






12 


BAHIA 


19/05/1 964 


6.40 






44.39 


58.70 


85.40 


109.02 


176.52 






13 


NABON 


22/03/1 966 


6.21 






44.03 


49.39 


76.51 


92.27 


174.11 






14 


PASTO CALLE 


06/10/1976 


5.92 






12.64 


17.82 


24.02 


43.76 


74.62 


103.96 




15 


ESMERALDAS 


09/04/1 976 


6.31 








15.96 


105.61 


173.95 


229.08 


278.52 




16 


TUMACO 


12/12/1 979 


7.70 




27.49 


68.94 


163.99 


250.60 










17 


BAEZA 


03/05/1 987 


7.00 




16.61 


28.70 


56.78 


80.47 


106.85 


132.32 


200.64 




18 


BAHIA DE 

carAquez 


08/04/1 998 


7.10 






350.53 


354.02 


369.14 


381.07 


422.51 


430.89 




19 


MACAS 


03/10/1995 


6.90 






160.93 


202.61 


222.90 


247.40 


257.29 






20 


PUJILf 


28/03/1996 


5.72 








33.57 


42.07 


60.69 


83.82 


118.03 


1 74.42 
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2. LEY DE ATENUACION DE INTENSIDADES 



Aguiar y Castro (2009) encontraron una ley de atenuacion de Intensidades al trabajar 
con la informacion macrosismica indicada en la tabla 3. Se trabajo con todos los sismos 
indicados en la tabla 3 sin clasificarlos en el sentido si estan asociados a fallamiento cortical o 
estan asociados al proceso de subduccion; clasificacion que se considera en el presente 
trabajo. En la tabla 3 la variable D es la distancia al foco. 

La ley de atenuacion de intensidades encontrada por Aguiar y Castro (2009) es la 
siguiente. 



7 = 3.3577 + 1.0013 -0.8856 Ln(D + 10) (1 ) 

Donde 7 es la intensidad en la escala de Mercalli Modificada; es la magnitud de 
Kanamori o momento sismico y 7) la distancia hipocentral. 



2.1 Reemplazo de curvas isosistas 



Ahora bien se volvio a estudiar el sismo de Ibarra del 16 de agosto de 1868, 
empezando por indicar que la magnitud es =7.7 y no 7.17 como consta en la tabla 3. 

Aguiar (2010). Este terremoto se denomina sismo de Ibarra, aunque en ciudades como 
Otavalo, Atuntaqui hubo mas dano que en Ibarra. Se trabajo ahora con la curva de isosistas 
del Observatorio Astronomico, indicada en la figura 2. 




Figura 2 Isosista del sismo de Ibarra de 1868. Fuente: Observatorio Astronomico. 

Tambien se cambio las curvas isosistas del sismo de 1906, por las proporcionadas por 
Gallego y Sarria ( 2006). La isosista de este sismo, con la que se trabajo en este articulo es la 
presentada en la figura 3. Este sismo es uno de los que mas energia a liberado a nivel mundial, 
en el Ecuador es el de mayor magnitud. 

Las curvas de isosistas del ultimo sismo que se cambio es la del sismo de Tumaco del 
12 de diciembre de 1979. La curva de isosista de este sismo es la indicada en la figura 4. 
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2.2 Nuevas curvas de isosistas 

Por otra parte se incorporo nuevas curvas de isosistas, que no fueron consideradas en 
ei estudio de Aguiar y Castro (2009). Las nuevas curvas corresponden a ios sismos de Manabi 
dei 13 de mayo de 1942 que causo mucho dano en edificios de hormigon armado ubicados en 
ia Ciudad de Guayaquii y se iiego a pensar de que ei hormigon armado no es un materiai 
adecuado para zonas sismicas pero con ei tiempo se vio todos ias deficiencias constructivas 
que se habian cometido, se encontro que en iugar de estribos ios hierros iongitudinaies de ias 
coiumnas habian side amarrados con sogas y cosas simiiares, de ahi que tuvieron un mai 
comportamiento. 





Figura 4 Isosista dei sismo de Tumaco de 1979. Puente: Gaiiego (2010) 



Revista CIENCIA Volumen 13 N.- 1 de 2010 



9 




Figura 5 Isosista del sismo de 1942. Fuente: Observatorio Astronomico. 

En base a la curva de isosistas se estimo que la profundidad focal del sismo de 1942 
fue de 50 km., la magnitud de este terremoto fue de 7.9. En la figura 5 se indica ias isosistas 
dei sismo de 1942. 




Figura 6 Sismo de Huila del 16 de novlembre de 1827. Fuente: Gallego (2010). 
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Se incorporo tambien al estudio el sismo de Huila del 16 de novlembre de 1827 que se 
registro en el Sur de Colombia y se lo sintio en el Ecuador. La Isosista de este evento se lo 
Indica en la figura 6. 

2.3 Sismos asociados a fallamiento cortical 

No solo que se camblaron las curvas de Isosistas y que se Incorporo nuevas curvas de 
otros terremotos, con relaclon al trabajo realizado por Aguiar y Castro (2009). SI no que, lo mas 
Importante es que se separaron los sismos de acuerdo a su origen. Se separo en sismos 
asociados a fallamiento cortical y sismos de subducclon. 

En la tabla 4 se presentan los 14 sismos asociados a fallamiento cortical con los cuales 
se obtuvo la sigulente ley de atenuaclon pero antes se destaca que son 69 puntos los que 
Intervinleron en el calculo de la ecuaclon. 

7 = 3.756 + 1.586 -1.861 * Ln(D + 10) (2) 

El significado de las variables ha sido ya Indicado. La desviaclon estandar de esta 
ecuaclon es 1 .21 . Garcfa (201 0). 



Tabla 4 Sismos asociados a fallamiento cortical. 



SISMOS CORTICALES 


# 


SISMOS 


FECHA 


Mw 


h 


D = (R''2 + h''2)''0,5 


X 


IX 


VIII 


VII 


VI 


V 


IV 


III 


II 


1 


RIOBAMBA 


04/02/1797 


7,70 


15,00 


25,26 


49,15 


88,69 


136,20 


182,40 


230,90 


276,77 


392,24 




2 


OTAVALO 


16/08/1868 


7,70 


33,00 


33,72 


40,53 


59,43 


68,54 












3 


PASTO CALLE 


15/09/1944 


5,92 


33,00 








34,03 


38,47 


45,15 








4 


PELILEO 


05/08/1948 


6,93 


10,00 


11,47 


14,80 


20,95 


28,94 


37,50 


53,90 


77,24 






5 


AMBATO 


05/08/1949 


6,37 


60,00 




77,79 


93,89 


143,76 




334,23 








6 


HUILA 


05/08/1949 


8,00 


33,00 




54,00 


105,54 














7 


NAPO 


11/05/1955 


6,98 


33,00 






40,06 


66,70 


180,72 










8 


PASA 


30/07/1960 


5,97 


21,00 








22,15 


24,29 


30,99 


55,86 


88,01 




9 


ALAUSI 


08/04/1961 


6,35 


24,00 








26,63 


43,15 


57,47 


83,58 






10 


NABON 


22/03/1966 


6,21 


42,00 






44,03 


49,39 


76,51 


92,27 


174,11 






11 


PASTO CALLE 


06/10/1976 


5,92 


7,40 






12,63 


17,82 


24,02 


43,76 


74,62 


103,96 




12 


BAEZA 


03/05/1987 


7,00 


10,00 




16,61 


28,70 


56,78 


80,47 


106,85 


132,32 


200,64 




13 


MACAS 


03/10/1995 


6,90 


160,20 






160,93 


202,61 


222,90 


247,40 


257,29 






14 


PUJILI 


28/03/1996 


5,72 


33,00 








33,56 


42,07 


60,69 


83,82 


118,04 


174,42 



2.3 Sismos asociados a ia subduccion 

En la tabla 5 se presentan los 8 sismos que fueron considerados para encontrar la ley 
de atenuaclon, asociados con la subduccion. Notese que se tienen 38 puntos. Al igual qup(§r] 
la tabla 4, en las ultimas columnas de la tabla 5 se indica la distancia al hipocentro de un 
cfrculo equivalente a la isosista obtenida. La Ley de Atenuaclon encontrada es la sigulente: 

7 = -1.484 + 1.402 -0.545 * Ln(D + 10) (3) 



Tabla 5 Sismos asociados a la subduccion. 



SISMOS SUBDUCCION 


# 


SISMOS 


FECHA 


Mw 


h 


D = (R''2 + h''2)''0,5 


X 


IX 


VIII 


VII 


VI 


V 


IV 


III 


II 


*r 


ESMERALDAS 


31/01/1906 


8,8 


25 


61,46 


81,06 


180,36 


212,94 


303,17 










*2 


MANABI 


13/05/1942 


7,8 


50 






128,97 


164,29 


217,67 










3 


BAHIA 


16/01/1956 


7,3 


25 






32,598 


53,006 


107,26 


214,95 


326,78 






4 


ESMERALDAS 


10/01/1958 


7,8 


60 








101,97 


146,88 


219,59 


284,28 






5 


BAHIA 


19/05/1964 


5,64 


34 






44,384 


58,699 


85,4 


109,02 


176,53 






6 


ESMERALDAS 


09/04/1976 


6,31 


10 








15,961 


105,61 


173,95 


229,08 


278,52 




*7' 


TUMACO 


12/12/1979 


7,7 


24 




28,37 


34,86 


89,05 


133,50 










7 


BAHIA DECARAQL 


08/04/1998 


7,1 


350,5 






350,53 


354,02 


369,14 


381,07 


422,51 


430,9 
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Las variables de la ecuacion (3) ya fueron indicadas; resta indicar que la desviacion 
estandar es 1 .16 



3. RELACION ENTRE ACELERACION MAXIMA E INTENSIDAD 



En las ecuaciones (2) y (3) se ha definido Leyes de Atenuacion de la Intensidad en 
funcion de la Magnitud y de la Distancia al hipocentro. Ahora, para fines de diseho sismo 
resistente se necesita definir una relacion con la cual se pase de Intensidad a Aceleraciones 
maximas del suelo. Esta relacion es de la siguiente forma: 

lOg^max (4) 

Donde es la aceleracion maxima esperada; es la intensidad en la escala de 

Mercalli Modificada; a y b son parametros que se obtienen en funcion de los dates con los 
cuales se trabajan. En la tabla 6 se indican los parametros a, b de algunos trabajos realizados. 



Tabla 6 Algunos trabajos sobre relacion entre Intensidades y Aceleraciones. 



Autor 


a 


b 


Trifunac y Brady (1975) 


0.300 


-0.014 


Bolt 


0.313 


0.340 


Murphy, O’ Brien (1977) 


0.250 


0.250 


Lomnitz (1974) 


0.333 


-0.500 


Saragoni et al. 1982 


0.345 


-0.222 


Saegesser y Baumbartner (1978) 


0.26 


0.19 


Medvedev-Sponheuer (1969) 


0.30 


-0.40 


0.30 


-0.1 



En las cinco primeras ecuaciones de la tabla 6, la Intensidad corresponds a 
(Mercalli Modificada), en las dos ultimas es (Medvedev, Sponheuer y Karnik). 

Estrictamente la ecuacion de Medvedev-Sponheuer (1969) es la siguiente: 

0.3 /„, - 0.4 < log < 0.30 /„„ -0.1 ( 5 ) 

Siete son los sismos con los cuales se hallo una relacion entre la aceleracion maxima y 
la Intensidad en la escala de Mercalli y son los registrados en Colombia en 1999, 2004 y 2008; 
en Chile en 2005, 2007 y 2010; en Peru en 2007. Las aceleraciones maximas del sismo de 
Chile de 2010 fueron proporcionadas por Boroschek et al. (2010) y son acelerogramas no 
corregidos. 

En el anexo 1 , se muestran los 99 datos de estos sismos con los cuales se hallo la 
siguiente relacion 



logUM«) = 0.3629 + 0.2355 »/„„ (6) 

Siendo la aceleracion maxima del registro; I la intensidad en la escala de 

Mercalli Modificada. En la figura 7 se muestran los pares de valores de intensidad y aceleracion 
maxima, tambien se encuentra graficada la ecuacion ( 6 ) hallada en el estudio. La desviacion 
estandar de esta ecuacion es 0.533 
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Intensidades [Imm] 



Figura 7 Relacion entre Aceleracion Maxima e Intensidad, hallada en el estudio. 



4. LEYES DE ATENUACION. PRIMERA PROPUESTA. 



Al reemplazar la ecuacion ( 6 ) en las ecuaciones ( 2 ) y ( 3 ) se tienen las nuevas 
Leyes de Atenuacion para el Ecuador, para sismos asociados a fallamiento cortical y para 
sismos de subduccion. Estas Leyes son las siguientes: 

• Fallamiento Cortical 

log(A/, ) = 1 .2474 + 0.3735 0.4383 *Ln{D + \0) ( 7 ) 

• Subduccion 

log(A^) = 0.0134 + 0.3302 -0.1284*Ln(D + 10) (8) 

Donde es la aceleracion horizontai en gais; ia magnitud y Z) ia distancia ai 

foco en km. La desviacion estandar de ias Leyes de Atenuacion encontradas, son: 0.65 para 
ios sismos asociados a faiiamiento corticai y 0.62 para ios sismos de subduccion. 

En ias figuras 8 y 9 se presentan estas ieyes para sismos de magnitud = 6 y 
Myy = 7, respectivamente. Tambien se presenta ia iey de atenuacion haiiada por Aguiar y 

Castro (2009). Dei anaiisis de estas figuras se desprende io siguiente: 1) La Ley de Atenuacion 
de Aguiar y Castro (2009) se aproxima bastante bien a ia Ley de Atenuacion encontrada en 
este estudio para sismos corticaies; 2) Es muy marcada ia diferencia entre ia Ley de 
Atenuacion para sismos corticaies y ia Ley de Atenuacion para sismos de subduccion. Esta 
uitima presenta vaiores muy bajos. 
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5. LEYES HALLADAS A PARTIR DE ACELEROGRAMAS 

Enia zona Central de Chile se han dado terremotos muy fuertes, con magnitudes del 
orden de 8.5, como los registrados en 1575, 1647, 1730, 1822, 1906, 1985 y 2010. Siendo el 
ultimo, el de mayor magnitud = 8.8 , mayor a la magnitud estimada del sismo Valdivia de 
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1906 que fue =8.5. Con esta pequena introduccion se pasa a describir brevemente el 

trabajo realizado por Ruiz y Saragoni (2005) para obtener leyes de atenuacion para sismos 
interplaca (zona de contacto de las placas tectonicas) e intraplaca (sismos asociados a 
fallamiento cortical). 

Ruiz y Saragoni (2005) encontraron las primeras Leyes de Atenuacion para Chile, para 
aceleracion, velocidad y desplazamiento, a partir de registros sismicos. Para el efecto 
trabajaron con 15 sismos con magnitud mayor a 6 {M ^ > 6) . En estos 15 eventos se tuvieron 

8 acelerogramas en roca dura >1500 (Suelo denominado SI) y 60 acelerogramas en 
suelo duro 1500 <V^ > 360 (Suelo S2). Donde es la velocidad de la onda de corte en m/s. 

Las ecuaciones encontradas para las aceleraciones horizontales (en gals), en un 

perfil de suelo S2, para fallamiento cortical es la ecuacion ( 9 ) y para sismos de subduccion 
son la (1 0) para un perfil de suelo SI y la (1 1 ) para un perfil S2. 



• Sismos intrapiaca (Faiiamiento Corticai) 



Ln{Af, ) = 8.2532 + 1 .2 * M - 2. 16 * Ln(D + 80) ( 9 ) 

• Sismos Thrust (Subduccion) 

) = 1.3863 + 1.3 -1.43*Ln(D + 30) (10) 

Ln(A/,) = 0.6931 + 1.28*M-1.09*Ln(D + 30) ) 



En la figura 10 se presentan las leyes de atenuacion para sismos de magnitud 6 y 7, 
encontradas aplicando las ecuaciones de Ruiz y Saragoni (2005) con las ecuaciones halladas 
en este estudio. Las graficas superiores corresponden a magnitud 6 y las inferiores a magnitud 
7; de igual manera las graficas de la izquierda son para sismos corticales y los de la derecha 
para sismos de subduccion. 

Para los sismos asociados a fallamiento cortical las leyes de atenuacion para Chile 
proporcionan valores superiores a las leyes de atenuacion encontradas en este estudio. Se 
destaca que la ley de atenuacion para Chile es para un perfil de suelo S2. 

Para los sismos asociados con la subduccion, la aproximacion entre las leyes de 
atenuacion de Ecuador y Chile tienen mejor aproximacion. Para distancias menores a 20 km., 
las ecuaciones de Ruiz y Saragoni (2005) proporcionan mayores valores; para distancias 
mayores a 20 km., sucede lo contrario la ley de atenuacion de Ecuador presenta mayores 
valores especialmente si se comprara con la Ley de Chile para un perfil de suelo SI . 
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Figura 10 Comparacion de leyes de atenuacion de Chile y Ecuador. Primera Propuesta 



6. LEYES DE ATENUACION PARA ECUADOR SEGUNDA PROPUESTA 



La segunda propuesta de Leyes de Atenuacion para Ecuador, consiste en utilizar la 
ecuacion propuesta por Saragoni et al. (1982) indicada en la tabla 6, que relaciona la 
aceleracion maxima con la Intensidad en la escala de Mercalli Modificada y que se la copia a 
continuacion. 



log!''™ ) = 0.345 0.222 (12) 

Se sustituye esta expresion, en las ecuaciones ( 2 ) y ( 3 ), encontradas en este estudio 
y se obtienen las Leyes de Atenuacion para el Ecuador de la segunda propuesta. Estas Leyes 
son las siguientes: 



• Fallamiento Cortical 

log [Af, ) = -1 .074 + QM1*M^~ 0.642 * ln(D + lO) (13) 

• Subduccion 

log ) = -0.734 + 0.484 * - 0. 188 * In (D + lO) (14) 



En la figura 11 se comparan las leyes de atenuacion halladas en este estudio, en la 
segunda propuesta, con las encontradas por Ruiz y Saragoni (2005). El formato de la figura 1 1 
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es similar al formato de la figura 10 per la que se omite la descripcion del contenido de las 
mismas. 





Figura 1 1 Comparacion de leyes de atenuacion de Chile y Ecuador. Segunda Propuesta. 



La Ley de Atenuacion hallada en la segunda propuesta se aproxima de mejor forma a 
la Ley hallada a partir de acelerogramas. Sin embargo las dos Leyes encontradas para el 
Ecuador son bastante buenas. 

7. CONCLUSIONES 

Lamentablemente en el Ecuador no se disponen de acelerogramas con valores 
mayores al 20% de la gravedad, con los cuales se podria encontrar una ley de atenuacion para 
ser utilizada en estudios de peligrosidad sismica empleando metodos deterministicos. Por esta 
razon se ha recurrido a informacion macrosismica que en forma general dan una buena idea de 
como se atenuo el movimiento del suelo con la distancia. 

Por otra parte, en base a 99 dates de aceleraciones maximas e intensidades de 7 
eventos sfsmicos registrados en los ultimos 12 ahos en: Colombia, Peru y Chile se hallo una 
relacion entre aceleracion maxima del suelo con intensidad, la misma que permitio pasar la Ley 
de Atenuacion de Intensidades a una Ley de Atenuacion de aceleraciones maximas del suelo. 
Con esta relacion se obtuvo una Ley de Atenuacion de la aceleracion horizontal del suelo para 
sismos asociados a fallamiento local y para sismos de subduccion 

Las ecuaciones encontradas son las indicadas en las ecuaciones (2) y (3). La ecuacion 
(2) encontrada en este estudio para fallamiento cortical se comparo con la ecuacion (1) que fue 
obtenida por Aguiar y Castro (2009) cuando no se clasifico la informacion macrosfsmica 
encontrado una muy buena aproximacion. Pero ahora tambien se tiene una Ley de Atenuacion 
para sismos de subduccion. 

Pero tambien se obtuvieron dos Leyes mas de Atenuacion, utilizando para el efecto la 
relacion de Saragoni et al. (1982) con la cual se paso de intensidad a aceleraciones 
horizontales. Estas dos nuevas Leyes de Atenuacion para el Ecuador, encontradas de 
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informacion macrosismica, fueron comparadas con leyes de atenuacion halladas en Chile pero 
a partir de acelerogramas, encontrando una muy buena aproximacion tanto para sismos 
asociados a fallamiento cortical como para sismos de subduccion. 

En resumen, en este trabajo se ha hallado una relacion que permite pasar de 
intensidades maximas a aceleraciones maximas, tema que se debe seguir trabajando en 
Raises que no tenemos suficientes acelerogramas de sismos fuertes, por este motive es que 
en el Anexo 1 se indica la base de dates. 

Se ha presentado dos propuestas para obtener Leyes de Atenuacion para el Ecuador, 
de las dos se recomienda la segunda ya que esta se aproxima mejor a las leyes de atenuacion 
halladas a partir de registros sismicos. Para esta segunda propuesta se tienen dos leyes una 
para fallamiento cortical y otra para subduccion. 
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ANEXO 1 

Datos con los cuales se hallo una relacion entre Aceleraciones e Intensidades 

^ Aceleraciones maximas(gals) 



Sismo 


Fecha 


Mw 


h 

(km) 


IX 


VIII 


VII 


VI-VII 


VI 


V-VI 


B 


IV 


III 


II 


Colombia 


15/11/2004 


7,2 


<33 






54,6 




















19,3 




















17,9 




















50,5 




















61,5 




















24,7 




















38,3 




















108,5 




















39,5 




















46,5 




















67,5 




















63,4 
















Chile 


13/06/2005 


7,9 


117,0 








259,0 


418,9 


66,7 


167,8 


66,0 


65,0 










211,9 


65,7 




183,4 


53,0 


128,0 










237,4 






157,9 


62,0 












709,3 






319,8 




















305,1 




















87,3 




















297,2 




















112,8 




















19,6 




















90,0 




















116,0 








Chile 


14/11/2007 


7,7 


60,0 




400,2 


91,0 




48,9 
















390,4 




82,8 




















74,1 




















59,8 












Peru 


15/08/2007 


7,9 


26 






272,0 




78,7 


58,4 














488,0 




54,1 


25,3 


















74,5 




















104,0 
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85,3 




















67,0 




















55,0 




















58,7 




















54,9 




















22,1 




















20,6 












Colombia 


24/05/2008 


5,9 


4,0 
















71,8 


3,5 


4,2 
















605,6 


3,4 


















48,2 




















37,7 




















6,5 




















23,9 




















6,4 




















4,6 




















4,3 




















11,8 




















7,1 






Chile 


27/02/2010 


8,8 


47,4 


462,1 


161,9 






90,4 














231,5 




















551,3 




















289,4 




















260,0 


















Colombia 


25/01/1999 


6,2 

(Ml) 


<15 




446,7 






8,0 






12,0 


1,9 






95,8 






6,8 








1,6 






73,6 














3,7 






99,1 














4,8 






63,3 














2,8 






97,4 














2,9 






28,0 














5,3 






25,6 














7,1 






57,5 




















182,2 



















Los datos del sismo de Colombia de 2004 fueron obtenidos de Salcedo y Alvarado 
(2007) y de Ojeda y Bermudez (2005). Para el sismo de Chile de 2005 las referencias son: 
Boroschek et al. (2006) y Salinas y Lazares (2005). Para el sismo de Chile de 2007, la fuente 
es Boroschek et al. (2007). En el sismo de Peru de 2007, la referenda es Tavera (2008). La 
fuente del sismo de Colombia de 2008 es INGEOMINAS (2008). La referenda del sismo de 
Chile de 2010 es Boroscheck et al. (2010). Finalmente, las citas bibliograficas del sismo que 
afecto al eje cafetero de Colombia en 1999 son: Gallego y Ordaz (1999); INGEOMINAS. 
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RESUMEN 



Se presentan dos modelos numericos de analisis sismico para puentes con aisladores 
de base tipo EPS (Frictional Pendulum System). En los dos modelos los pilotes son modelados 
con elementos finitos lineales con masa uniforme distribuida; los estribos con elementos finitos 
rectangulares concentrando las masas en puntos discretes seleccionados; el suelo que rodea a 
los pilotes es modelado con resortes con comportamiento no lineal en base a la curva carga 
deformacion del suelo. 

En el modelo 1 no se considera el la presencia del suelo que se encuentra en la parte 
exterior de los estribos y solo se considera que el EPS trabaja como un resorte de rigidez . 
En el segundo modelo se toma en cuenta el suelo que esta en los lados de los estribos, 
mediante una malla de elementos finitos y los EPS son modelados como resortes de rigidez 

y amortiguamiento viscose . La rigidez k^ fue obtenida en funcion del perfodo objetivo de la 
estructura, del radio de curvature del EPS y de la carga que gravita sobre el EPS. El factor de 
amortiguamiento fue obtenido experimentalmente de las curves de histeresis de un EPS. 

Los dos modelos fueron comparados realizando el analisis sfsmico del Puente Sur, que 
forma parte del grupo de puentes construidos sobre el estuario del rfo Esmeraldas con 
aisladores de base EPS. 

Por otra parte se presentan los nuevos programas que ban side implementados para el 
software CEINCI-LAB. 



ABSTRACT 



Two numeric models of seismic analysis are presented for bridges with insulators of 
base type EPS (Frictional Pendulum System). In the two models the piles are modeled with 
lineal finite elements with distributed uniform mass; the stirrups with rectangular finite elements 
concentrating the masses on selected discreet points; the floor that surrounds to the piles is 
modeled with springs with non-linear behavior based on the curve it loads deformation of the 
floor. 
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In the pattern 1 are not considered the presence of the floor that is in the external part 

k 

of the stirrups and alone it is considered that the FPS works as a spring of rigidity. * In the 
second model he/she takes into account the floor that is in the sides of the stirrups, by means of 

k 

a mesh of finite elements and the FPS they are modeled as springs of rigidity * and viscous 
c k 

reduction. * The rigidity was * obtained in function of the objective period of the structure, of 
the curvative radius of the FPS and of the load that gravitates on the FPS. The factor of 

reduction was * obtained experimentally of the curves of hysteresis of a FPS. 

The two models were compared carrying out the seismic analysis of the South Bridge 
that is part of the group of built bridges on the estuary of the river Emeraldas with base 
insulators FPS. 

On the other hand the new programs are presented that have been implemented for the 

CEINCI-LAB software. 



1. DESCRIPCION DE PUENTE 



El Puente Sur, es otra de las estructuras que el Cuerpo de Ingenieros del Ejercito ha 
construido en el estuario del rfo Esmeraldas con aisladores de base tipo pendulo de friccion 
FPS (Frictional Pendulum System). Es un puente que tiene dos estribos tipo cajon en sus 
extremes y dos pilas centrales, cubriendo una luz de 160 m.; en la figura 1 se muestra una vista 
longitudinal del puente, se aprecia que las luces exteriores son de 50 m., y la luz central es de 
60 m. 



VISTA LONGITUDINAL DEL PUENTE SUR 

BAJA PENDENTE = 2,10 %; LONG. = 160.00 m 



6.30 1.70 



1 .70 3.30 




ESTRIBO DERECHO 



ESTRIBO IZQUIERDO 



Figura 1 Vista longitudinal del puente Sur. 



El estribo izquierdo se halla sobre una cimentacion de hormigon ciclopeo de 2.0 m., de 
alto, con una dimension en planta de 12.0 m., por 21.0 m., aproximadamente ya que la forma 
es trapezoidal como se ve en la figura 1. Sobre el hormigon ciclopeo se encuentra una losa 
maciza de hormigon armado de 0.70 m., de alto, que vendrfa tambien a ser parte de la 
cimentacion; las dimensiones de esta viga en planta son de 8.20 m., por 19.10 (redondeando 
las cantidades). Sobre esta losa maciza viene el cajon propiamente dicho cuyas paredes 
interiores son de 0.30 m., y las paredes exteriores son de 0.35 m. En el sentido de la dimension 
larga existen 3 paredes interiores y dos paredes exteriores; en el sentido de la dimension corta 
las paredes son de 0.30 m. Estas paredes son las que forman las celdas del estribo que luego 
son rellenadas con material petreo. En la figura 2, a la izquierda se muestran las paredes cajon 
a las que se hace referenda y a la derecha se muestra la viga trapezoidal sobre la que se 
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apoyan los aisladores FPS, para el modelo numerico de calculo la viga trapezoidal se ha 
convertido en una vIga rectangular. 

El estrlbo derecho es similar al estrlbo Izquierdo pero con dimenslones mas pequehas, 
asi la losa que se encuentra sobre la base de hormigon ciclopeo es de 0.60 m., de alto. La 
altura total del estrlbo derecho es 10.15 m., mientras que la altura total del estrlbo Izquierdo 
esta por los 13.0 m., en los dos casos se ha Incluldo los 2 m., de alto de hormigon ciclopeo que 
esta bajo tierra. 

Los dos apoyos centrales que se observan en la figura 1, tienen la misma concepcIon 
estructural pero con dimenslones diferentes. Cada uno de estos apoyos descansa sobre 10 
pllotes pre-barrenados de 1.22 m., de diametro, hincados a 17 m., de profundidad y colocados 
en forma circular en un diametro de 14 m. 



Sobre los pllotes se encuentra la denomlnada viga cabezal, que en realldad es una 
losa maciza de 1 .50 m., de alto, y tiene un diametro de 14 m. Sobre esta viga cabezai se halla 
la plla circular que es maciza en una altura de 1 .0 m., tiene un diametro de 6 m., luego esta plla 
se convierte en una plla hueca con un diametro Interior de 4.80 m., de tal manera que las 
parades de la plla tienen un espesor de 0.60 m.; en la parte superior esta plla es trapezoidal y 
al final se tiene una viga superior de las sigulentes dimenslones en planta: 3.00 X 16.60 m., el 
alto de esta viga es de 1 .0 m. 



O 
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Figura 2 Vista en planta y un corte interior de estrlbo izquierdo. 



En la figura 3 se describe en forma grafica el puente Sur, en la parte superior izquierda 
se aprecia uno de los estribos tipo cajon, durante su construccion, en el cual se han colocado 
los tres aisladores FPS. En la parte superior derecha se aprecian las parades del estrlbo, en su 
parte posterior. 

En las fotograffas de la figura central se observe las 6 vigas de acero tipo I que 
soportan la superestructura del puente. Los FPS estan en la mitad de estas vigas 
longitudinales, claro esta que sobre los FPS se encuentra una viga transversal tipo I que se 
denomina diafragma. En la fotograffa central derecha se ve mejor que el FPS se encuentra en 
la mitad de las vigas longitudinales. 
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En la parte inferior de la figura 3 se observan las pilas centrales del puente, a la 
izquierda en el proceso constructive y a la derecha cuando estaba terminada. A la izquierda se 
ven las seis vigas de acero que cubren dos luces exteriores de 50 m., y una luz central de 60 
m. En la fotografia inferior derecha se aprecia el sitio en el cual se unieron las vigas para cubrir 
la luz de 160 m., se aprecia en la viga de acero la junta vertical. 

En la figura inferior derecha se aprecia los cabezales de los pilotes que tienen una 
forma circular, la pila circular de 6 m., de diametro exterior y la forma como se reduce esta 
seccion hasta llegar a la viga superior de 3.0 m por 1 6.60 m. 

Antes de presenter aspectos del modelo numerico de calculo, es importante presenter 
aspectos constructivos que se deben tener en cuenta en estructuras con aisladores de base. 
Lo importante es tener en cuenta que los aisladores se van a mover bastante (del orden de los 
20 a los 30 cm., ante un sismo muy fuerte) y que no deben encontrar impedimentos durante 
este movimiento. Se recuerda que la disipacion de energfa durante un sismo severo se va a 
dar en los EPS. Aguiar ef a/2008., Constantinou (1990), Almazan (2001). 





Figura 3 Descripcion grafica del puente Sur. 
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A la izquierda de la figura 4 se ve la viga de acero, una de las parades del estrlbo cajon 
y en la parte superior se encuentra la losa del tablero del puente. Con una flecha se senala que 
entre la pared del estrlbo y la losa existe una separaclon de mas de 3 cm., esto es para que 
pueda moverse llbremente el tablero y no exista fricclon con las paredes del cajon del estrlbo. 

Ahora, a la derecha de la figura 4 se ve que entre el tablero del puente y el estrlbo se 
ha dejado una junta de construcclon de 40 cm., dimension Inusual con respecto a los puentes 
que se construyen sin alsladores de base. Nuevamente esto es debldo a la presencia del 
alslador, cuyo funclonamlento se Indica en el proximo apartado. La tapa junta es una placa de 
acero. 




2. FUNCIONAMIENTO DE LOS FPS 



Existen tres generaclones de FPS, el primero data de los ahos noventa en el cual se 
tiene una superficle concava, un desllzador y una superficle recta. Zayas et al. (1989). La 
segunda generaclon de los FPS es de finales del siglo XX y de principles del siglo XXI, en el 
cual se tiene doble curvature y un desllzador. Fenz y Constantinou (2006). 

Los de la tercera generaclon, corresponde al utlllzado en los puentes de Esmeraldas, 
son dos alsladores de doble curvature en uno solo dispositive. En la figura 5 se aprecia uno de 
ellos, se ven las superficies concaves que son de acero Inoxidable muy blen pulldas para 
permitir del desllzamlento del dispositive Interior que tamblen tiene dos superficies de acero 
concavas y en la parte central se tiene el desllzador. Ante movimlentos sfsmicos de pequeha 
Intensidad trabaja unicamente el alslador Interior y ante movimlentos sfsmicos muy fuertes 
trabajan los dos desllzadores el Interior y el Interior de este. Los alsladores vienen recublertos 
por una goma que Impide la entrada de polvo que podrfan Impedir el movimlento del desllzador. 




Figura 5 Alslador FPS utlllzado en los puentes de Esmeraldas 

En la figura 6 se pretende expllcar el funclonamlento de los FPS, en sentido transversal 
al eje del puente, por ello se ha dibujado un estrlbo con los tres FPS. Se presentan dos flechas, 
la superior corresponde a las fuerzas generadas por el movimlento sfsmico y la Inferior son las 




26 



Revista CIENCIA Volumen 13 N.- 1 de 2010 



fuerzas que se generan en la superficie de contacto del desllzador. En otras palabras cuando 
hay un sismo el desllzador se mueve pero Internamente se generan fuerzas de sentido 
contrarlo que se oponen al movimlento. En forma muy elemental se puede declr que estas 
fuerzas de fricclon valen v N . Donde V es el coeficlente de fricclon que existe entre el 
desllzador y la superficie concava de acero; este coeficlente es funclon de la velocldad de 
deformaclon y de la preslon. N es la fuerza normal; para pequehas deformaclones se tiene 
que la fuerza normal es Igual al peso que gravlta sobre el EPS. 



En este artfculo se modela el FPS como un resorte horizontal que tiene una rigidez 



la misma que se halla de la sigulente manera: I) se determine la carga vertical que gravlta en el 
alslador; II) se encuentra el radio de curvature del FPS, con la sigulente expreslon: 



T = 27T 




(1 ) 



Donde T es el perfodo objetivo que se desea tenga el sistema de alslaclon; R es el 
radio de curvatura; g es la gravedad. Ill) se halla la rIgIdez con la formula sigulente. Aguiar 
etal. (2008). 



kb 



m g 
R 



( 2 ) 



Slendo m la masa que se obtlene dividlendo el Peso que se transmite al alslador para 
la gravedad g . 



.18^90, 




Figura 6 Esquema de trabajo, en sentido transversal, de los FPS de un estribo. 



El perfodo objetivo T se impone el calculista. En el case de Esmeraldas este perfodo 
fue de 3 s.; Se trabaja con perfodos que estan entre los 2 y 3 segundos, ya que para esa zona 
del espectro de diseho sfsmico las ordenadas espectrales son bajas. 

Al ser bajas las ordenadas espectrales; la fuerza que se transmite a la estructura por 
efecto del sismo es baja; esta es una de las ventajas de utilizer aisladores de base, en general. 
Otra de las ventajas es que al trabajar con perfodos altos se esta flexibilizando a la estructura 
con lo cual el alslador filtra las altas frecuencias del sismo, las que producen daho. 
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3. NUEVOS PROGRAMAS DE CEINCI-LAB 



Con el desarrollo de esta investigacion se han incorporado nuevos programas a 
CEINCI-LAB que son especfficos para el anallsis sismico del puente Sur de Esmeraldas pero 
que pueden ser llgeramente modificados y complementados para el anallsis estatico o sismIco 
de otro tlpo de estructuras, por este motive se Indica su uso de algunos de ellos. 

3.1 Generacion de mallas en vigas continuas 



Lo mas compllcado en elementos finitos es la generacion de la matriz que contlene a 
los vectores de colocaclon de una estructura. Por este motive se desarrollo un programa 
denomlnado vc_acero_tres_vanos_ef que obtlene las coordenadas generallzadas de cada 
uno de los nudos y la matriz de vectores de colocaclon de cada uno de los elementos para una 
viga continua de 3 vanos. 

En este programa los nudos se numeran de Izquierda a derecha y de abajo hacia 
arriba, lo propio sucede con los elementos. Por otra parte los apoyos son rodlllos, de tal 
manera que existe un grade de llbertad horizontal en cada apoyo. 

Como esta orlentado al anallsis sismico, se concentre las masas en los ejes de los 
apoyos como se aprecia en la parte Inferior de la figura 7. De tal manera que las coordenadas 
principales son las horizontales y que estan asocladas con las masas concentradas. 
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COORDENADAS GENERALES 
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COORDENADAS PRINCIPAEES 



Figura 7 Grades de llbertad del programa vc_acero_tres_vanos_ef 

Una vez que se numeran las coordenadas principales, se pasa a numerar todos las 
demas coordenadas horizontales de los nudos, como se Indica en la parte superior de la figura 
7, luego se numeran las coordenadas verticales. 

El programa vc_acero_tres_vanos_ef reporta: la matriz que contlene a los vectores 
de colocaclon de cada uno de los elementos; el numero total de grades de llbertad y las 
coordenadas generallzadas de cada uno de los nudos. 

El usuarlo de CEINCI-LAB solo tiene que IndIcar el numero de elementos fInItos en 
direcclon X y en direcclon Y, de un solo vano. El programa asume que este mismo numero de 
elementos se consideran en los otros vanos. En la figura 7 se han considerado 4 elementos 
finitos en direcclon X y 3 elementos finitos en direcclon Y. Por clerto las dimenslones Indicadas 
de los vanos corresponden a la vIga del Puente Sur; la vIga es de acero de 2.47 m., de alto. 



2.47 
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3.2 Matriz de rigidez y de masas de las vigas de acero 



El programa que halla la matriz de rIgIdez y de masas de las vigas de acero del puente 
Sur se denomlna rigidez_viga_acero_PS. El programa primero halla la matriz de rigidez de la 
viga continua, con todos los grados de llbertad; luego condense la matriz de rigidez a las 
coordenadas laterales Indicadas en la parte Inferior de la figure 7. Aguiar (2004). La matriz de 
masas encuentra tamblen para el mismo sistema de coordenadas principales. El llstado del 
programa es el sigulente. 

function [KA,MA]=rigidez_viga_acero_PS 

% 

% Programa para encontrar la matriz de rigidez y de masa para anaiisis 
% sismico de viga superior de acero dei Puente Sur. 

% 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%% 

% XX XX XX XX 

% 

% 50 m. ! 60 m. ! 50 m. ! 

% 

% Roberto Aguiar 

% CEINCi-ESPE 

% Febrero de 201 0 

% 

% Matrices de Rigidez y de Masas de Viga Superior de Acero 

% 

E=20000000.0; % Modulo de elasticidad del acero en T/m2 

[X,Y,XP,YP]=coord_ef(1 80,3,1 60.0,2.47); %Coordenadas de Elementos Finitos 

[VC,ngl,cg]=vc_acero_tres_vanos_ef(60,3); 

for i=1 :180*3; espe(i)=0.045*6; end % espesor medio por 6 vigas 

[K]=rigidez_estructura_ef(VC,180,3,XP,YP,E,espe,0.20,ngl); 

% Condensacion de matriz de rigidez 

for i=1 :16; r1(i)=i; end; for i=1 :ngi-16; r2(i)=16+i; end 

Kaa=K(r1 ,r1 );Kab=K(r1 ,r2);Kbb=K(r2,r2); 

KA=Kaa-Kab*inv(Kbb)*Kab'; 

MA=zeros(16,16); %Matriz de masas 

masa=1 60*2.47*0.045*7.85/9.8 ;MA(1 ,1 )=M A(1 ,1 )+masa/1 6; 

MA(2,2)=MA(1,1);MA(3,3)=MA(1,1);MA(4,4)=MA(1,1); 

MA(5,5)=MA(1,1);MA(6,6)=MA(1,1);MA(7,7)=MA(1,1);MA(8,8)=MA(1,1); 

M A(9,9)=MA(1 ,1 );M A(1 0,1 0)=MA(1 ,1 );M A(1 1,11 )=MA(1 ,1 ) ;M A(1 2,1 2)=MA(1 ,1 ); 

M A(1 3,1 3)=MA(1 ,1 );M A(1 4,1 4)=M A(1 ,1 );MA(1 5,1 5)=M A(1 ,1 );M A(1 6,1 6)=MA(1 ,1 ); 
return 
% — end— 



Como era de esperarse el programa llama a otros programas de la llbrerfa de CEINCI- 

LAB. 

3.3 Matriz de rigidez y de masas de estribo izquierdo 



En la figura 8 se IndIcan las 28 coordenadas principales considerados para el estribo 
Izquierdo. El programa rigidez_estribo_izquierdo_PSurDETALLADO encuentra la matriz de 
rigidez y de masas. Se trabajo con 13 Elementos FInItos (EF) en X y 2 EF en Y, en la 
cimentaclon de hormigon ciclopeo; el modulo de elasticidad empleado para este hormigon fue 

E = 1264911 T I . Para la losa de hormigon armado que se halla encima y que tiene una 
altura de 0.70 m., se calculo con 9 EF. 
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Continuando, para arriba, se tiene el estribo formado solo por las paredes de hormigon 
armado, en una altura de 4.95 m., y un ancho de 7.28 m. Para este sector se considero 8 EF 
en sentido X y 5 EF en sentido Y, en total 40 EF. 

Para la viga que esta bajo el EPS se calculo con 2 EF en X y 2 EF en Y, en total 4 EF. 
En la parte lateral Izquierda de esta vIga existen paredes de hormigon armado, este sector fue 
modelado con 6 EF en X y 2 EF en Y. 

Luego hay todo un sector que tiene solo paredes de hormigon armado, esto fue 
calculado con 6 EF en X y 3 EF en Y. Las dimenslones de este sector del estribo es de 6.40 
por 2.83 m. FInalmente se tiene la losa superior del estribo que constituye el tablero del puente 
y tiene 0.25 m., de alto, en este sector se considero 6 EF en X y 1 en Y. El modulo de 

elasticldad del hormigon fue E = 1870828 T I . 




Figura 8 Coordenadas principales de estribo izquierdo. 



3.4 Matriz de rigidez y de masas de apoyos 

Como se habfa manifestado hay dos apoyos muy semejantes pero diferentes en sus 
dimenslones, para cada uno de ellos se elaboro un programa de computacion para encontrar la 
matriz de rigidez y de masas; rigidez_apoyo1_central_PS es el nombre del primer programa 
para el apoyo que esta mas proximo al estribo izquierdo y rigidez_apoyo2_central_PS es el 
otro programa. El numero de coordenadas principales del apoyo 1 es 56 y del apoyo 2 es 52. 
En los dos casos la ultima coordenada considerada como principal es una vertical asociada con 
la componente de desplazamiento vertical del punto en el cual se asienta el EPS. Esto con el 
proposito de a future incorporar modelos mas sofisticados de analisis. 



3.5 Matriz de rigidez y de masas de estribo derecho 



El programa rigidez_estribo_derecho_PS encuentra la matriz de rigidez y de masas 
en coordenadas principales del estribo derecho. En este case se tienen 22 coordenadas 
principales, las mismas que se indican en la figura 1 0. 
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Para los apoyos centrales, los pilotes fueron modelados con EF lineales con masa 
concentrada y para los estrlbos con EF rectangulares con dos grados de llbertad por nudo. El 
marco teorico se presenta en el articulo anterior. Aguiar y Garcia (2010). 
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Figura 10 Coordenadas principales de estribo derecho. 



4. ANALISIS SISMICO LONGITUDINAL CON MODELO 1 



En la figura 1 1 se presentan los grados de libertad con los cuales se realiza el analisis 
sismico en sentido longitudinal. Se tienen 176 coordenadas principales entre las que se 
incluyen dos coordenadas verticales en los apoyos centrales y son las coordenadas 84 y 136. 

En este articulo el EPS es modelado como un resorte de rigidez , hallado con la 

ecuacion ( 2 ) y que actua en los siguientes grados de libertad 159, 160, 161 y 162. En un 
future trabajo se modelara al EPS como un elemento pequeno, que esta entre los apoyos y la 
viga. 



El suelo que rodea a los pilotes es modelado, con resortes que actuan en los grados de 
libertad 29 a 68 para el apoyo 1 y 85 a 124 para el apoyo 2. La curva carga deformacion 
utilizada para el suelo es la indicada en la figura 12. Aguiar et a! (2010), Beltran (2008). Se 
entiende que a partir de los 7 m., la curva carga deformacion es la misma. 

Con CEINCI-LAB se indica solo las propiedades del suelo de un eje de pilotes y el 
programa gener a g dl resortes obtiene las propiedades de los demas ejes de pilotes de un 
apoyo. De tal forma que se debe dar como datos dos ejes de pilotes, el uno correspondiente al 
apoyo 1 y el otro eje de pilotes corresponds al apoyo 2. 

La estructura en la forma analizada tiene un periodo fundamental de 1.34 s. Antes de 
finalizar la explicacion del modelo 1 y pensando en el modelo 2 que se indica posteriormente se 
debe manifestar lo siguiente. El modelo 1 no considera la presencia del suelo que se encuentra 
al lado de los estribos en el movimiento sfsmico que si considera el modelo 2. La otra diferencia 
es que el modelo 1 no considera el amortiguamiento del EPS que si considera el modelo 2. Se 
destaca una vez mas que se esta presentando un modelo aproximado para el comportamiento 
del EPS. 
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Figura 1 1 Grades de libertad considerados en el analisis sfsmico longitudinal 




El sismo de analisis se indica en ia figura 13 y corresponde a un aceierograma artificiai 
que tiene una duracion de 20 s., y que es compatibie con ei espectro dei CEC-2000 para un 
perfii de sueios S3. A ia derecha de ia figura 13 se muestra ei espectro de respuesta eiastico 
haiiado para 5 % de amortiguamiento. 
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Figura 13 Acelerograma artificial y su espectro elastico. 



5. RESULTADOS CON MODELO 1 



En la figura 14 se muestra la respuesta en desplazamlentos, en el estrlbo Izquierdo; 
apoyo 1; apoyo 2; estrlbo derecho. En cada una de estas subestructuras se selecclonaron 4 
puntos los mismos que estan Identificados por los grades de llbertad (GDL). De acuerdo al GDL 
en la figura 1 1 se puede ver la ubicaclon del punto selecclonado. 






15 t(s.) 




Estrlbo Izquierdo 



Apoyo 1 






15 t(s.) 




Apoyo 2 



Estribo Derecho 



Figura 14 Respuesta de desplazamlentos hallados con modelo 1 . 



Se presenta la respuesta a nivel de cimentaclon, en los estrlbos o a nivel de los pllotes 
en los apoyos centrales; ese es un punto. Luego Interesa conocer los desplazamlentos bajo los 
alsladores y sobre los alsladores, finalmente es Interesante conocer los desplazamlentos en la 
parte superior de la viga de acero. En la figura 14 se aprecia que los desplazamlentos en la 
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cimentacion o en los pilotes, son extremadamente bajos; lo propio sucede con los 
desplazamientos laterales bajo el FPS. Sobre los FPS los desplazamlentos son altos y en la 
parte superior de la viga se mantlenen estos desplazamlentos; es declr, los desplazamlentos 
relativos del tablero con respecto al alslador son practicamente nulos. 

La respuesta a nivel de fuerzas horizontales se IndIcan en la figura 15, para los mismos 
grados de llbertad de la figura 14. Los cortantes sobre el alslador son bajos, el valor maximo es 
del orden de las 10 T., (era preocupaclon encontrar el cortante en este punto para ver la 
demanda sismica que llega a los pernos de anclaje superior de los FPS). En la parte Inferior de 
los FPS los cortantes son mas altos pero ahf se tiene una gran vIga de hormigon armado (los 
pernos de anclaje de los FPS estan Inmersos en esta gran viga). 
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Figura 15 Respuesta de Fuerzas con Modelo 1 



15 «(s.) 



En la figura 16 se muestra la respuesta maxima en valor absoluto, en la cara externa 
del estribo Izquierdo (graficos de la parte superior) y en la cara externa del primer apoyo 
(graficos de la parte Inferior). Se presenta la respuesta de desplazamlentos, a la Izquierda y de 
fuerzas a la derecha. 



Tanto los estrlbos como las pllas practicamente no se desplazan horizontalmente hasta 
donde se encuentran los FPS, son los alsladores los que se desplazan notablemente y luego 
nuevamente aparece una recta vertical en los graficos que estan a la Izquierda de la figura 16, 
lo que significa que no hay desplazamlento. Este comportamlento es debldo a que con la 
presencia de los FPS el peso de la superestructura detlene los desplazamlentos de la 
subestructura. SI no habrfa FPS los desplazamlentos de las pllas centrales serfa considerable. 
Curves simllares se obtuvieron para el apoyo 2 y el estribo derecho. 
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Figura 16 Respuestas maximas en estribo izquierdo y apoyo 1 . 



1 Adovo 1 


^ 







X: 75.49 
Y:18 



10 20 30 40 50 60 70 F(T.: 



6. ANALISIS SISMICO LONGITUDINAL CON MODELO 2 



En la figura 17 se muestran los dos modelos simplificados, utilizados para los EPS, a la 
izquierda se indica el modelo 1, descrito anteriormente y a la derecha el modelo 2, en que se 
incorpora el amortiguamiento del EPS, como viscose. 



Kb 



Figura 17 Modelos simplificados de los aisladores FPS 

Cb =2 <^b ^ ) 

Donde es el factor de amortiguamiento del aislador; es la masa que gravita 
sobre el aislador y kb la rigidez del FPS. 

En la figura 18 se presenta la curva de histeresis de un aislador FPS, en el eje de las X 
se presentan los desplazamientos y en el eje de las Y, la fuerza aplicada. El factor de 
amortiguamiento se obtiene de esta curva con la siguiente ecuacion. 





4 



Atu E, 



(4) 



Siendo la energfa disipada por el sistema inelastico y es igual area de un cicio de 
histeresis; E^ es la energfa elastica que absorbe el sistema. Aguiar (2003). En la figura 18, 
proporcionada por Landazuri (2009) se observe que el factor = 0.29 , para el aislador 
ensayado. Para ser un poco conservador se trabajo con = 0.25 . 
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Figura 18 Curva de histeresis de un aislador FPS. 

El modelo 2, incorpora tambien el suelo que se encuentra en la parte exterior de los 
estribos. Para esto se trabajo con un EF rectangular con dos grados de libertad por nudo, 
cuyas funciones de forma son las indicadas en la tabla 1 . 



Tabla 1 Funciones de forma para elemento finite rectangular. 



Coordenada natural 


Deformada 

elemental 


Funcion de Forma 


1 


ul 


(Z)i =0.25(-5 + l)(-r + l) 


vl 


(Z), =0.25(-5 + l)(-r + l) 


2 


u2 


(z)2 = 0.25(5 + l)(-r + l) 


v2 


^2 = 0.25(5 +iX-f+i) 


3 


u3 


(^3 = 0.25(5 -Fl)(r + l) 


v3 


<2)3 = 0.25(5 -Fl)(r + l) 


4 


u4 


<2)4 =0.25(-5 + l)(r + l) 


v4 


<2)4 =0.25(-5 + l)(r + l) 



Donde s,t son las coordenadas naturales del EF que varfa entre -1 y 1. Notese que 
no se ha incorporado en el modelo la presion del agua, para ir paso a paso mejorando y 
ensehando los EF.Cook et al. (2002). 



7. NUEVOS PROGRAMAS DE CEINCI-LAB 

En la figura 19 se muestra el modelo numerico de calculo del estribo izquierdo, 
incorporando la presencia del suelo. Las masas concentradas corresponden a las coordenadas 
principales indicadas en la figura 8. A la izquierda de estas masas se ha colocado el suelo con 
sus respectivos grados de libertad. 
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Figura 19 Coordenadas generales consideradas en el modelo del estrlbo Izquierdo. 



Para generar la matriz que contiene a los vectores de colocacion del suelo se ha 
generado el programa vc_suelo_ef. Los dates de entrada de este programa son: el numero de 
elementos finitos en sentido X; el numero de elementos finitos en sentido Y; una matriz que 
contiene las coordenadas generalizadas del primer eje de vertical de la malla y una matriz que 
contiene las coordenadas generalizadas del ultimo eje vertical de la malla; en esta matriz 
primero se coloca el grade de libertad horizontal y luego el vertical, en una fila. 

Por cierto vale la pena recordar que el vector de colocacion se forma con los grades de 
libertad del elemento empezando desde el vertice inferior izquierdo y se continua en forma anti 
horaria, VC = [ 0 0 0 0 261 262 259 260], para el primer elemento de la figura 19. La 
primera fila de la malla indicada en la figura 19 no es factible generar con vc_suelo_ef. 

Las coordenadas principales del estrlbo izquierdo se indica en la figura 20. Notese que 
se ha considerado como coordenada principal el desplazamiento vertical bajo el aislador, por 
esto se tiene ahora 29 grades de libertad. El programa que encuentra las matrices de rigidez y 
de masa del estrlbo izquierdo es: rigidez_estribo_izquierdo_suelo_PSurDETALLADO . En 
este programa se considero que el modulo de poisson del suelo es 0.4 y que el modulo de 
elasticidad del suelo se halla con la siguiente ecuacion. 

_ £„(i-k-2k") 



l-v 



( 4 ) 
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Donde es el modulo de deformaclon del suelo, se trabajo con =300T I ^ 
V es el coeficlente de polsson del suelo. 




El programa rigidez_estribo_derecho_suelo_PS obtlene la matriz de rigidez y de 
masa en coordenadas principales del estrlbo derecho. En la figura 21 se IndIcan las 
coordenadas totales del estrlbo Incluyendo el suelo a la derecha y en la figura 22 las 
coordenadas principales. En este caso no hay como apllcar directamente el programa 
vc_suelo_ef, para encontrar los vectores de colocaclon del suelo, pero se puede IndIcar el 
vector de colocaclon de los elementos que estan en contacto con el estrlbo derecho y los 
demas grados de llbertad se pueden generar con vc_suelo_ef. 

Notese que tamblen se ha considerado como coordenada principal, el desplazamlento 
vertical en el apoyo. 



8. COORDENADAS PRINCIPALES DE MODELO 2 



En la figura 23 se IndIcan los grados de llbertad con los que se realizan el anallsis 
sfsmico considerando el modelo 2. Se tiene en total 176 grados de llbertad, como coordenadas 
principales. El perfodo de vibraclon de la estructura es de 1 .34 s. 

La masa vertical provenlente de la superestructura se concentre en los grados de 
llbertad 29 (estrlbo Izqulerdo), 85 (apoyo 1), 137 (apoyo 2) y 160 (estrlbo derecho). Estos son 
grados de llbertad verticales que fueron considerados como coordenadas principales. En cada 
uno de ellos la masa se obtuvo dividlendo el peso (170 T * 3) para la gravedad. La carga 
vertical en cada EPS es 170 T. 
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Figura 21 Coordenadas totales del estribo derecho con el suelo. 
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Figura 22 Coordenadas principales on estribo derecho con Incorporaclon del suelo. 
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Figura 23 Grades de libertad de modelo con interaccion suelo estructura en estribos. 



9. RESULTADOS CON MODELO 2 



Se hace notar al lector que el sistema de coordenadas del modelo 1 , es diferente al 
sistema de coordenadas del modelo 2. Para el nuevo modelo la historla de desplazamlentos se 
Indica en la figura 24 para los grades de libertad Indicados que son: a nivel de cimentaclon, 
bajo el FPS, sobre el FPS y en la parte superior de la viga de acero. 
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Figura 24 Respuestas de desplazamiento con Modelo 2. 
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Las respuestas de desplazamiento son mucho menores a las obtenidas con el modelo 
1, se reducen a menos de la mitad como se podra apreclar al comparar la figura 24 con la 
figura 14. No coinclden los grados de llbertad pero los puntos que se estan graficando son los 
mismos en los dos modelos. 




Apoyo 2 Estribo derecho 

Figura 25 Respuestas en el tiempo de las fuerzas horizontales. Modelo 2. 
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Figura 26 Respuestas maximas encontradas con el Modelo 2. 
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En la figura 25 se muestra la historla de las fuerzas halladas en el estrlbo Izquierdo, 
apoyo 1 , apoyo 2 y estrlbo derecho. Son diferentes a las encontradas en el modelo 1 y que 
constan en la figura 1 . 

Por ultimo en la figura 26 se presentan los valores maximos en las cara exterior del 
estrlbo Izquierdo y en la cara exterior del apoyo 1 . En este grafico es evidente la diferencia que 
existe entre estos resultados y los del modelo 1 que se encuentran en la figura 26. Notese que 
en el modelo 2 los desplazamlentos maximos estan por el orden de los 20 cm., en camblo en el 
modelo 1 estos desplazamlentos estan por los 50 cm. 



10.CONCLUSIONES 



Se ha presentado dos modelos simpllficados para el anallsis sfsmico de puentes con 
alsladores de base EPS. En el primer modelo no se considero la presencia del suelo que se 
encuentra en la parte exterior de los estrlbos, en el segundo si se lo considero medlante un 
modelo de elementos finitos rectangulares. 

En el primer modelo el EPS es modelado como un resorte de rigidez ; en el segundo 
modelo a mas del resorte se considera el amortiguamlento del alslador c^, el factor de 
amortiguamlento del alslador fue obtenido experlmentalmente. Para comparar los dos 

modelos se realize el anallsis sfsmico del puente Sur, que se encuentra sobre el estuarlo del rfo 
Esmeraldas, es un puente de 160 m., de luz conformado por dos estrlbos y dos apoyos 
centrales. 

Para cada uno de los modelos Indicados se han desarrollado programas muy 
especfficos que se Incrementan a la llbrerfa del software educative y practice CEINCI-LAB, por 
este motive es que se han descrito con clerto detalle ya que pueden utlllzarse en otros puentes 
simllares con pocos camblos. 

La principal conclusion del estudio es que no Incorporar el amortiguamlento del EPS en 
el modelo numerico de calculo conduce a obtener desplazamlentos laterales que son 
aproximadamente el doble de los que se hallan en el modelo 2 que considera el 
amortiguamlento de los EPS. Lo propio sucede a nivel de fuerzas. 

Con el desarrollo de este artfculo, en forma didactica, se aspira haber aportado al 
desarrollo de la Ingenlerfa SIsmo Resistente. 
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RESUMEN 

La industria de los pesticides naturales en la actualidad a tornado una gran 
importancia, es asi que en algunos paises se ha tornado la decision de reemplazar a los 
pesticides sinteticos por los naturales, El uso indiscriminado de pesticides sinteticos ha 
acarreado una serie de problemas en la salud humane, asi como tambien han repercutido 
gravemente en la flora y la fauna animal, ya que estos pesticides no son biodegradables y con 
el consume continue y prolongado causan series enfermedades letales a los organismos 
vivientes. En nuestro pafs se esta empezando a desarrollar tecnologfas para la fabricacion de 
pesticides naturales. 

En este artfculo se trabajo en la obtencion de los pesticides naturales en base a los 
principios actives del Eucalipto y del Ajo, que llevandolos a soluciones homogeneas en 
diversos medios de solucion, se llegaron a establecer las condiciones mfnimas y optimas para 
el control y la mitigacion de plages en algunas plantas ornamentales de la ESPE asf como 
tambien en cultivos de limon y tomate de arbol, es asf que en forma separada o en forma 
conjunta utilizandose ya sea el eucalipto o el ajo, se logro mitigar al Planococcus citri, 
Conchuela y al Pulgon. 



ABSTRACT 

The natural pesticide industry at present taken on great importance, so that in some 
countries have taken the decision to replace synthetic pesticides by natural Indiscriminate use 
of synthetic pesticides has led to a series of problems in the human health and also have 
severely affected the flora and fauna animal, because these pesticides are not biodegradable 
and continuous long-term consumption cause serious diseases lethal to living organisms. Our 
country is beginning to develop technologies for the production of natural pesticides. 

This article is working on obtaining natural pesticides based on the active principles of 
Eucalyptus and Garlic, leading to homogeneous solutions in various means of settlement, it 
came to establish minimum and optimum conditions for the control and mitigation pest in some 
ornamental plants in the ESPE as well as lemon and tomato crops tree, so that separately or 
together using either eucalyptus or garlic, it was possible to mitigate Planococcus citri, 
Conchuela and Aphid . 



1. INTRODUCCION 

El eucalipto tiene muchas propiedades beneficas para la salud humana y animal asf 
como tambien para el control de plagas debido a su composicion es decir sus principios 
activos. El principal components del aceite esencial es el eter oxido terpenico cineol o 
eucaliptol, constituyendo el 70-80%, ademas contiene hidrocarburos monoterpenicos (10% de 
b-pineno, g-pineno, p-cimeno, camfeno, a-felandreno, b-felandreno, limoneno, mirceno), 
hidrocarburos sesquiterpenicos (1 % de aromadendreno, b-cariofileno, a-gurjuneno), alcoholes 
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monoterpenicos (borneol, linalol, terpineno, a-terpineol, 1,2% de transpinocarveol, geraniol), 
alcohol sesquiterpenico (6,9% de globulol; epiglobulol, trazas de eudesmol, ledol, virldiflorol), 
aldehidos allfaticos (butirico, valerlanico, caprolco), citronelal, carvona, acetate de citronelllo, 
acetate de geranllo, 3,6% de acetate de a-terpinll (15-20). 

El ajo y sus componentes, come son los sulfoxides: S-metll-L-cIsteina, S-propenll-S- 
clsteina, S-glutatl6n, g-glutamll-S-alllcIsteina, g-glutamll-S-alll-mercapto-L-cIsteina, ademas se 
encuentran sales mlnerales (selenlo), azucares, lipidos, amlnoacldos esenclales, saponosidos, 
terpenos, enzimas, flavonoldes y otros compuestos fenollcos, contlene acelte volatll (1-3%) 
sulfurado (disulfuro de alllo), sulfoxide de S-alll (L)-(-i-)-clstema o allfna (0,2%), vltamlnas A y C, 
nicotllamida (7,8). La allfna se transforma en allcina (ver figura 1) y otros compuestos 
azufrados (tiosulfinatos), por la acclon de la enzima alllnasa. Estos ultimos son muy Inestables 
y se transforman con extrema rapidez en otros compuestos organosulfurados: sulfuro de dialllo, 
disulfuro de dlalllo(mayorltarlo en la esencia de ajo),trlsulfuro dedlalllo y ajoenos, todos ellos 
solubles en medio oleoso. 



Se considera que 1 mg de allfna equivale a 0,45mg de allcina. El ajo posee 
propledades antloxidantes, tiene efectos antiterogenico, antihepatotoxico y anticancerfgeno, 
propledades Inhibitorlas de la agregaclon plaquetarla, la allcina es un potente agente 
antibacterlano frente a numerosos organismos patogenos gram-positivos y gram-negativos se a 
verlficado su actividad sobre bacterlas de los generos Staphylococcus, Streptocuccus, Vibrio y 
Bacillus (12,15). Los aceltes esenclales ejercen una potente acclon antlseptica, balsamica y 
expectorante en la vfa respiratorla, antifungica, antituberculosis, antiviral (2,8,11,12,14), 
antiheimfntico, antlespasmodico, vasodllatador(4), diuretico, hipcolesterolemico, estimula la 
secreclon estomacal y blllar. 

En este estudio se prepararon los pesticldas naturales en base a los principlos activos 
del Eucallpto y el Ajo, los cuales fueron llevados a diferentes soluclones homogeneas y a 
concentraclones adecuadas para determlnar la ellmlnaclon de plagas como son el Planococcus 
citrl, Conchuela y el Pulgon. En diferente tlpo de plantas ornamentales de la ESPE entre las 
cuales tenemos la Schefflera, Mlllonarla, Palmera Washington, Mala Madre, Helecho y en 
cultivos de LImon y Tomate de arbol. 

2. METODOLOGIA 

2.1 Preparacion de los principlos activos de las plantas de Eucallpto y Ajo 

En el procedimlento de preparacion se recolectaron las hojas de la planta de eucallpto 
tenlendose culdado de que estas fueran lo mas frescas posibles y llbres de parasites para 
poder someterlas rapidamente al proceso de extracclon (18,19). Para esto se recolectaron 4 
sacos llenos de hojas de eucallpto, se procedio a fracclonarlas en pedazos mas pequehos, 
seguldamente se las coloco en un reciplente de acero ahadlendose agua de manera que las 
hojas ocuparan la mitad del volumen de agua, seguldamente se tape el reciplente el cual 
estaba provisto en su parte superior de una cabeza clalsen unida a un refrigerante de cobre, se 
calento a temperatura de ebulllclon del agua, y se recoglo por destllaclon 200 ml de acelte puro 
de eucallpto. 




Aliina 



Alicina 



Figura 1 Transformacion enzimatica de la aliina en alicina. 
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Para la extraccion de los principios actives del aceite esencial del ajo, se aisle el agente 
active basice del aje, la allina (6-8) que cuande es liberada interactua cen una enzima llamada 
allinasa y de esta ferma se genera la alicina, la sustancia que centiene el eler caracteristice y 
penetrante del aje. Para la ebtencion de la alicina se trabaje a temperatura y presion ambiente, 
para le cual en primer lugar se descascare el aje, seguidamente se le semetio a un precese de 
trituracion, luege de licuacion, para luege extraer el principle active cen alcehel etilice al 96%, 
en extraccienes sucesivas de hasta 5 veces. Para ebtener el principle active mas pure se 
cencentre el alcehel etilice ebteniendese una masa gelatinesa de un celer amarille parde, entre 
etras sustancias tenemes: sulfoxides S-metil-L-cisteina y S-propenil-S-cisteina, S-glutation, g- 
glutamil-S-alil cisteina, g-glutamil-S-alil-mercapto-L-cisteina. 



2.2 Elaboracion de soluciones como vehiculos para los principio activos 

Una vez obtenidos los principios activos tanto del eucalipto como del ajo, se procedio a 
la elaboracion de productos en soluciones homogeneas. Se probaron varias formulaciones 
entre las cuales las mas efectivas fueron: 

F1 : La cual esta constituida per el principio active del aceite de eucalipto en diferentes 
proporciones, con el vehiculo en el cual tenemos etanol y propilenglicol. F2: Constituida por el 
principio active del aceite de eucalipto en diferentes proporciones, con el vehiculo en el cual 
tenemos etanol, propilenglicol y lanolina hidrosoluble. F3: Formada por el principio active del 
aceite de eucalipto en diferentes proporciones, con el vehiculo constituido por etanol, 
propilenglicol, lanolina y emulsion de silicona al 60% con el objeto de que el topico se fije mas 
permanente en la planta sea en epoca seca o Iluviosa. F4: Teniendo por principio active del ajo 
en diferentes proporciones, con el vehiculo formado por etanol, propilenglicol y lanolina 
hidrosoluble. F5: Constituida por el principio active del aceite de ajo en diferentes proporciones, 
con el vehiculo formado por etanol, propilenglicol, lanolina hidrosoluble y emulsion de silicona al 
60% para una mayor fijacion. F6: Formada por la combinacion del principio active del aceite de 
eucalipto y el principio active del ajo ambos en diferentes proporciones, en vehiculo formado 
por etanol, propilenglicol, lanolina hidrosoluble y emulsion de silicona al 60%. F7: La cual tiene 
en su composicion el extracto concentrado del principio active del ajo esto es una solucion 
amarillenta muy espesa obtenida por la condensacion de varias semanas de la fase final en la 
extraccion con alcohol etflico al 96%, se utiliza como vehiculo el etanol, propilenglicol, lanolina 
hidrosoluble y emulsion de silicona al 60%. 



2.3 Aplicacion sobre sistemas biologicos 

En esta fase se determine la aplicacion de los productos en la concentracion mfnima 
necesaria en base a los principios activos del eucalipto o del ajo, en diferentes plantas 
ornamentales de la ESPE entre las cuales tenemos, la Schefflera, Millonaria, Palmera 
Washington, Mala Madre, Helecho y en cultivos de Limon y Tomate de arbol, para la mitigacion 
de las plagas o sistemas biologicos encontradas las cuales fueron el Planococcus citri, 
Conchuela y el Pulgon. 



2.4 Analisis de ef-.-ectividad de los principios activos sobre los sistemas 
biologicos 

Para observar la efectividad de la aplicacion de los principios activos del eucalipto y el 
ajo llevados en diferentes soluciones homogeneas, sobre plagas se utilize un estereoscopio 
(14,15) marca LEICA S4E provisto con placa de asentamiento para cajas petri asf como 
tambien se utilizo el analisis en situ (29,30) sobre las plantas ornamentales de la ESPE y en los 
cultivos de limon y tomate de arbol infestadas con plagas. 
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3. RESULTADOS Y DISCUSION 

3.1 Preparacion y aplicacion de los principios actives del Eucalipto para la 
eliminacion de plagas en la planta Schefflera 

Para la aplicacion sobre sistemas biologicos primeramente se tomb el aceite de 
eucalipto 22,34g con vehiculo de etanol y propilenglicol para llegar a obtener la cantidad 
necesaria y suficiente de 1000,00g, esta formulacion se aplico a la especie de planta Schefflera 
la cual se encuentra atacada por la plaga Planococcus Citri, Cotonet o tambien llamada 
Cochinilla algodonosa(27-30). Esta plaga parte de la savia que toma la excreta como liquido 
azucarado brillante (melaza) sobre el cual se asienta el hongo la negrilla. Esta clase de 
ectoparasites pertenecen a la clase de aracnidos son rastreros atacan a las plantas ricas en 
savia en la cual se hospedan hasta destruirla completamente y secarla en un ataque masivo y 
sin control, el planococcus citri puede destruir a toda la planta en un periodo comprendido de 
varies semanas, ademas de acuerdo a las observaciones en el laboratorio de esta plaga se 
determine que pueden desplazarse a grandes distancias buscando su medio de desarrollo, por 
lo que estos tambien se encuentran en la tierra aledaha a las plantas, la hembra deposita los 
huevos de forma ovalada y color amarillo en un saco ovigeno algodonoso de forma irregular. 
Los machos se originan de nifas al completar su segundo estadio, presenta de 2 a 4 
generaciones anuales que se superponen por lo que se puede encontrar en cualquier epoca 
del aho en cualquiera de sus estados. Sin embargo en el invierno disminuye su densidad y es 
mas probable encontrar masas de huevos. Como se puede observer en la figure 2 al 
planococcus citri c ' ' ' .... 



A esta concentracion, actuo sobre la colonia del Planococcus citri, pero no en toda la 
poblacion. Para probar la efectividad del producto se recolectaron muestras de hojas y tallos 
atacados por el planococcus citri, una parte se coloco en cajas petri y otra parte se las dejo a la 
intemperie, ambas rociadas con la solucion FI que igualmente ataco a una parte de la 
poblacion. 

Se aumento la concentracion de los principios activos del eucalipto a 41 .99g por litro de 
vehiculo, se cojieron 20ml de esta solucion y se aforo a II, como resultados tenemos que se 
elimino paulatinamente la colonia de Planococcus citri, es decir atacandose a su sistema 
operative produciendose una paralizacion del animal seguido de un proceso de secamiento. 

Con la misma concentracion anterior de los principios activos del eucalipto a 41.99g 
por litro de vehiculo, se cojieron esta vez 40ml de esta solucion y se aforo a 1 1, como resultados 
tenemos que se elimino rapidamente a la colonia de Planococcus citri, es decir atacandose a 
su sistema operative produciendose una paralizacion del animal seguido de un proceso de 
secamiento, como se puede observer en la figure 3, siendo esta una aplicacion muy efectiva 




Figura 2 Planococcus Citri atacando una hoja de Schefflera 
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para el control del ectoparasite tanto por apllcaclon por aspersion dIrecta en la planta asi como 
tamblen en las pruebas efectuadas en el laboratorlo en caja petrl observadas al estereoscoplo. 




Figura 3 Proceso de secamiento en la eliminacion del planococcus citri 



La efectividad de esta soluclon se probo en varlas plantas scheffleras ubicadas en 
diferentes dependencies de la ESPE, en las cuales tenemos el edificio central, en el 1, 2, y 3 
piso, asi como tamblen en el edIfIcIo nuevo en el primer piso. 



3.2 Aplicacion del pesticida obtenido del Eucalipto en otras plantas 



La concentraclon minima requerlda obtenida en base a los principles actives del 
Eucalipto de 41 ,99g por litre de vehiculo, en la cual se cojieron 20ml de esta soluclon y se aforo 
a II, fue apllcada a varlas plantas, en las cuales se observe el nivel de mitigaclon del 
Planococcus citrl y tamblen de las conchuelas. En la tabla 1 se puede observer la efectividad 
del pesticide obtenido en base a los principles actives del Eucalipto para el control de la plage. 



Tabla 1 Pesticida en base de eucalipto a la concentraclon de 41 .99g/l y 20ml a 1 1, sobre 
planococcus citri 



Tipo de planta 


Principios activos del 
Eucalipto 

Grado de concentraclon 


Tipo de plaga 
PLANOCOCCUS CITRI 
Resultado 


Schefflera 


41.99g/ly20mla1l 


Repelida 


Millonaria 


41.99g/ly20mla1l 


Repelida 


Palmera Washington 


41.99g/ly20mla 11 


Repelida 


Mala Madre 


41.99g/ly20mla1l 


Repelida 


Helecho 


41.99g/ly20mla1l 


Repelida 


Limon 


41.99g/ly20mla1l 


Repelida 


Tomate de arbol 


41.99g/ly20mla1l 


Repelida 



Como se puede observer en la tabla 1 la accion de los principios actives del eucalipto 
sobre las diferentes plantas es mas bien repelente, es por eso que se cojieron 40ml de la 
soluclon anterior como se puede ver en la tabla 2. 
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Tabla 2 Pesticida en base de eucalipto a la concentracion de 41.99g/l y 40ml a 11, sobre 
planococcus citrl 



Tipo de pianta 


Principios activos del 
Eucalipto 

Grado de concentracion 


Tipo de piaga 
PLANOCOCCUS CITRI 
Resultado 


Schefflera 


41 ,99g/l V 40ml a 11 


Eliminacion rapida 


Mlllonarla 


41.99g/l y 40ml all 


Eliminacion rapida 


Palmera Washington 


41 ,99g/l V 40ml a 11 


Eliminacion rapida 


Mala Madre 


41 ,99g/l y 40ml a 11 


Eliminacion rapida 


Helecho 


41 ,99g/l V 40ml a 11 


Eliminacion rapida 


LImon 


41 ,99g/l V 40ml a 11 


Eliminacion rapida 


Tomate de arbol 


41 ,99g/l V 40ml a 11 


Eliminacion rapida 



La concentracion minima efectiva con accion mitigante sobre ia coionia de pianococcus 
citri se ia puede apreciar en ia tabia 2. 



3.3 Preparacion y aplicacion de los principios actives del Ajo en la 
eliminacion de plagas en diferentes plantas 

En ia apiicacion dei pesticida obtenido a partir de ios principios activos dei ajo, se 
probaron varias concentraciones, de ias cuaies ia mas efectivas son ias que se prepararon a 
partir de 370g de ia soiucion etanoiica dei ajo con sus principios activos, ademas utiiizando 
como vehicuio propiiengiicoi, ianoiina anhidra y soiucion de siiicona ai 60% en un totai de 493g, 
de este preparado se separaron 50 mi y se aforo a un iitro, como se puede apreciar en ia tabia 
3. 



Tabia 3 Pesticida en base de ajo a la concentracion de 370g/493g y 50ml a II, sobre 
planococcus citrl 



Tipo de pianta 


Principios activos del Ajo 
Grado de concentracion 


Tipo de piaga 
PLANOCOCCUS CITRI 
Resultado 


Schefflera 


370g/493g y 50ml a 1 1 


Repellda 


Mlllonarla 


370g/493gy50mla1l 


Repellda 


Palmera Washington 


370g/493g y 50ml a 1 1 


Repellda 


Mala Madre 


370g/493g y 50ml a 1 1 


Repellda 


Helecho 


370g/493g y 50ml a 1 1 


Repellda 


LImon 


370g/493g y 50ml a 1 1 


Repellda 


Tomate de arbol 


370g/493g y 50ml a 1 1 


Repellda 



En la tabia 3 se observe que a la concentracion indicada ei pianococcus citri no es 
eiiminado sino soiamente repeiido por io que es muy probabie de que en ia pianta se produzca 
una nueva reinfestacion, entonces es necesario aumentar ia concentracion efectiva como se 
indica en ia tabia 4, ia cuai eiimina rapidamente ai pianococcus citri. 



Tabia 4 Pesticida en base de ajo a la concentracion de 370g/493g y 100ml a II, sobre 
planococcus citrl 



Tipo de pianta 


Principios activos del Ajo 
Grado de concentracion 


Tipo de piaga 
PLANOCOCCUS CITRI 
Resultado 


Schefflera 


370g/493gy lOOmlall 


Eliminacion rapida 


Mlllonarla 


370g/493g y 1 00ml a 1 1 


Eliminacion rapida 


Palmera Washington 


370g/493gy lOOmlall 


Eliminacion rapida 


Mala Madre 


370g/493g y 1 00ml a 1 1 


Eliminacion rapida 


Helecho 


370g/493g y 1 00ml a 1 1 


Eliminacion rapida 


LImon 


370g/493g y 1 00ml a 1 1 


Eliminacion rapida 


Tomate de arbol 


370g/493g y 1 00ml a 1 1 


Eliminacion rapida 
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3.4 Analisis de la efectividad de los principios actives concentrados del ajo 

Seguidamente se prepare una solucion mas concentrada del pesticida del ajo para ello 
se utilize, 557g de la solucion etanolica del ajo, propllengllcol, etanol, lanollna, emulsion de 
slllcona y se puso 1.5g del concentrado pure de la solucion etanolica del ajo en 17g de agua, 
para completar los 681 .5g del total de la solucion, la cual fue apileada con buenos resultados, 
como se Indica en la tabla 5. 

Tabla 5 Pesticida en base de ajo a la concentracion de 557g y 1.5g/681.5g sin afore, sobre 



planococcus citri 



Tipo de planta 


Principios activos del Ajo, 
solucion etanolica y 
concentrada, 

Grado de concentracion 


Tipo de plaga 
PLANOCOCCUS CITRI 
Resultado 


Schefllera 


557g V 1 .5g/681 .5g sin aforo 


Eliminacion rapida 


Millonaria 


557g V 1 .5g/681 .5g sin aforo 


Eliminacion rapida 


Palmera Washington 


557g y 1 .5g/681 .5g sin aforo 


Eliminacion rapida 


Mala Madre 


557g V 1 .5g/681 .5g sin aforo 


Eliminacion rapida 


Helecho 


557g V 1 .5g/681 .5g sin aforo 


Eliminacion rapida 


Limon 


557g V 1 .5g/681 .5g sin aforo 


Eliminacion rapida 


Tomate de arbol 


557g y 1 .5g/681 .5g sin aforo 


Eliminacion rapida 



3.5 Aplicacion del pesticida formado per los principios activos del ajo y del 
eucalipto 

Se preparo el pesticida formado por el eucalipto y el ajo, para lo cual en la formula de 
370g solucion etanolica de ajo, se anadieron 8g de principios activos del eucalipto, ademas 
utilizando como vehiculo propllengllcol, lanollna anhidra y solucion de slllcona al 60% en un 
total de 501 g, para combatir al planococcus citri, como se puede observer en la tabla 6. 



Tabla 6 Pesticida en base de ajo y del eucalipto a la concentracion de 370g y 8g/501g sin aforo, 
sobre planococcus citri 



Tipo de planta 


Principios activos del Ajo, y 
del Eucalipto 
Grado de concentracion 


Tipo de plaga 
PLANOCOCCUS CITRI 
Resultado 


Schefflera 


370g y 8g/501g sin aforo 


Eliminacion rapida 


Millonaria 


370g V 8g/501g sin aforo 


Eliminacion rapida 


Palmera Washington 


370g V 8g/501g sin aforo 


Eliminacion rapida 


Mala Madre 


370g V 8g/501g sin aforo 


Eliminacion rapida 


Helecho 


370g V 8g/501g sin aforo 


Eliminacion rapida 


Limon 


370g y 8g/501g sin aforo 


Eliminacion rapida 


Tomate de arbol 


370g V 8g/501g sin aforo 


Eliminacion rapida 



Como se puede ver en la tabla 6, la efectividad del producto de la combinacion de los 
principios activos tanto del eucalipto como del ajo es de una eliminacion rapida del planococcus 
citri. 



3.6 Mitigacion del pulgon en base a los principios activos del ajo 

En las algunas plantas cheflera, tomate de arbol, y millonaria ademas del planococcus 
citri se encontro el pulgon, como se puede apreciar en las figuras 4 y 5, a los cuales se los 
conoce tambien como afidos (21), se caracterizan por su rapida multiplicacion y cuando atacan 
una planta puede verselos cubriendo en su totalidad las hojas y brotes, son insectos 
chupadores, provistos de un largo pico mediante el cual absorben la savia de las plantas y se 
presentan en variedad halada y sin alas. Los pulgones son los principales transmisores del 
virus, pican en una planta infectada y al picar en otra sana le inyectan el virus. 
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Figura 4 Puigon verde Figura 5 Pulgon amarillo 



Las variedades del pulgon son resistentes a las concentraclones Indicadas 
anterlormente en las diversas pruebas efectuadas, ya que solamente los repelia pero no los 
ellmlnaba, por lo que se probo solamente con el pesticlda del ajo, como se puede observer en 
la table 7, formado por los principlos activos en soluclon etanolica, el cual ellmino rapidamente 
el pulgon, y ademas quedaba una capa aceltosa adherlda a las superficies en contacto de las 
plantas la cual permitia que por algunos dias no se acerque el pulgon. 



Tabla 7 Pesticlda a base de ajo en soluclon etanolica 



Tipo de planta 


Principlos activos del Ajo en 
soluclon etanolica 
Grade de concentracion 


Tipo de plaga 
PULGON 
Resultado 


Schefflera 


500ml puro 


Eliminacion rapida 


Millonaria 


500ml puro 


Eliminacion rapida 


Tomate de arbol 


500ml puro 


Eliminacion rapida 



4. CONCLUSIONES 

Las plantas de eucalipto y ajo ban demostrado ser unos excelentes pesticides naturales, esto 
es en base a los resultados obtenidos, ya que estos pesticides no atacan solamente a un tipo 
de plaga sino que tienen un amplio margen de accion plaguicida, debido a sus principlos 
activos, los cuales actuan en el organismo del animal paralizando sus funciones vitales, pero 
no atacando a la planta en ninguna de sus funciones, dentro de las concentraclones mfnimas 
requeridas. 

Del eucalipto el principio active con accion plaguicida a demostrado ser el 1,8-Cineol, para el 
ajo el principio active con mayor accion plaguicida es la Allina o sulfoxido de S-alil-cistefna 
(aminoacido azufrado), que se transforma en Alicina por accion de la enzima Allinasa, estos 
principlos activos actuando en forma aislada o conjunta son muy buenos plaguicidas. 

Con los pesticidas obtenidos del ajo y del eucalipto, se probaron en varies plantas entre las 
cuales tenemos: Schefflera, Millonaria, Palmera Washington, Mala Madre, Helecho. Limon, 
Tomate de arbol, las mejor concentracion de accion plaguicida para el eucalipto es de 41 ,99g/l 
y 40ml a II, para el ajo la mejor concentracion de accion plaguicida es de 370g/493g y 100ml a 
II, para la soluclon etanolica del ajo y su extracto concentrado es de 557g y 1.5g/681.5g sin 
aforo, cuando se utilizan para accion plaguicida los principlos activos del ajo y del eucalipto 
juntos su mejor concentracion encontrada es de 370g y 8g/501g sin aforo y solamente los 
principlos activos puros del ajo en soluclon etanolica. 

En las plantas Schefflera, Millonaria, Palmera Washington, Mala Madre, Helecho. Limon, 
Tomate de arbol, se encontraron principalmente las plagas Planococcus Citri y el Pulgon, para 
el case del Planococcus citri se observe que este era vulnerable tanto para el eucalipto como 
para el ajo en las concentraclones de efectividad indicadas, pero en el case del pulgon este 
presentaba gran resistencia tanto para el eucalipto como para el ajo, e incluso si se trabajaba 
con los dos pesticidas juntos, por lo que de acuerdo a los ensayos realizados el unico que 
actuo en el sistema nervioso del animal paralizando sus funciones biologicas fue el extracto de 
los principlos activos del ajo en soluclon etanolica, ya que este a mas de atacar al animal 
formaba una capa aceltosa adherents y protectora del ataque del pulgon en la planta. 
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Una vez realizada la curacion de plagas en las plantas descritas, con cualquiera de los 
pesticldas menclonados ya sea del eucallpto o del ajo o de los dos juntos, es necesarlo repetir 
por lo menos cada quince dias la apllcaclon para evltar una reinfestaclon ya que los huevos de 
estos ectoparasites se pueden encontrar en los sectores aledanos a la planta sobre todo en la 
tierra. 



5. RECOMENDACIONES 



Seguir Investigando el uso benefico de los pesticldas naturales obtenidos de diversas plantas 
en la mitigaclon de plagas, para no depender del uso de los pesticldas sinteticos que son 
toxicos para el hombre y el medio amblente. 
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RESUMEN 

Se evaluo la eficiencia de reduccion de nitrogeno amoniacal, fosforo total y demanda quimica 
de oxigeno (DQO), en agua residual sintetica y semi-sintetica utilizando una cepa axenica de 
Nostoc sp. y un consorcio microbiano, inmovilizados en perlas de alginato de calcio y en 
suspension. Tras un periodo de 21 dias de tratamiento, se observe, para el agua residual semi- 
sintetica que el cultivo inmovilizado de Nostoc sp. presento la mayor eficiencia de reduccion de 
nitrogeno y fosforo, mientras que para el agua residual sintetica, no existio una reduccion de 
nitrogeno con respecto a la concentracion inicial del mismo con ninguno de los tratamientos, y 
la mayor reduccion de fosforo fue observada en el cultivo inmovilizado del consorcio 
microbiano. Transcurridos 18 dfas, se observe que para la DQO, en ambos ensayos, la mayor 
eficiencia de reduccion la presentaron los cultivos en suspension, tanto del consorcio como de 
la cepa axenica, incluyendo tambien el tratamiento sin inoculo, para el caso del agua residual 
semi-sintetica (que solo conto con la microbiota propia del agua residual). Finalmente se 
realizaron pruebas de toxicidad en Daphnia pulex con el agua tratada correspondiente al cultivo 
de cianobacterias axenico inmovilizado, sin evidenciarse mortalidad en los organismos de 
prueba. 



ABSTRACT 

It was evaluated the efficiency of ammonia, total phosphorus and chemical oxygen demand 
(COD) reduction in synthetic and semi-synthetic wastewater using an axenic strain of Nostoc 
sp. and a microbial consortium immobilized in calcium alginate beads and in suspension. After a 
period of 21 days of treatment, it was observed for the semi-synthetic wastewater that the 
immobilized culture of Nostoc sp. had the highest reduction efficiency of ammonia and 
phosphorus. For synthetic wastewater, there wasn’t ammonia reduction on the concentration 
with none of the initial treatments, and further reduction of phosphorus was observed in the 
cultivation immobilized microbial consortium. Within 18 days, it was observed that for CQD, in 
both trials, semi-synthetic and synthetic wastewater, that the greater efficiency of reduction was 
introduced by the suspension cultures, both the consortium as axenic strain, also including the 
control treatment, in the case of semi-synthetic wastewater (which only had the original 
microbial charge of wastewater). Finally toxicity tests were carried out on Daphnia pulex, using 
the treated water of the cultivation of axenic cyanobacteria immobilized without mortality evident 
in organisms test. 
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1. INTRODUCCION 

Las cianobacterias, o algas verde-azules, son un grupo diverse y amplio de organismos 
procariotes Gram-negativos que pertenecen al Dominio Bacteria, presentan una gran 
importancia debido a sus capacidades fotosinteticas y de fijacion de nitrogeno. Poseen una 
amplia distribucion ecologica, ocupando un rango de habitats que inciuyen desde zonas bastas 
oceanicas hasta desiertos o aguas termaies (Herrero et al., 2001; Haii et al., 1982). Son 
organismos fotoautotroficos que asimiian CO 2 a traves dei Cicio de Caivin en presencia de iuz, 
aunque aigunas especies son capaces de reaiizar crecimiento quimioheterotrofico en ia 
compieta oscuridad 0 a bajas intensidades de iuz (Herrero, etal., 2001 ; Fay, 1992). 

Ei agua residuai representa un medio apropiado para ei crecimiento de cianobacterias 
(Abaide et al., 1995), pues posee bicarbonatos, nitrates y fosfatos, ios cuaies son tres aniones 
para ios que estas han desarroiiado mecanismos aitamente eficientes para tomar y retener 
(Ernst et al., 2006) razon por ia que han side consideradas como una aiternativa para ei 
tratamiento de aguas residuaies (Chacon etal., 2004; Abaide et al., 1995). 

Ei cuitivo de cianobacterias en aguas residuaies ofrece dos ventajas simuitaneas, ia 
primera, ei tratamiento de dichas aguas y, por otra parte, ia posibiiidad de producir biomasa 
que posteriormente podria ser expiotada en forma de compiementos proteinicos, aditivos 
aiimenticios 0 farmaceuticos, energia como biogas y combustibies, fertiiizantes entre otros 
vaiiosos quimicos (Maiiick, 2002). 

Se ha comprobado que su use en ei tratamiento de aguas residuaies es potenciaimente 
provechoso en ia reduccion de contaminantes como ia DQO, DBO5, soiidos suspendidos 
totaies, metaies pesados, entre otros (Ei-Bestawy, 2008). Asi tambien ei diseho de sistemas de 
cianobacterias inmoviiizadas, representa una opcion de tratamiento economicamente eficiente 
(Maiiick & Rai, 1993). 

En uitimos estudios, se ha evidenciado su efectividad tambien en ia remocion de 
metaies traza dei agua, ampiiando asi ios campos de apiicacion de ias cianobacterias en 
materia de biorremediacion de aguas (Baptista & Vasconceios, 2006). 



2. MATERIALES Y METODOS 

2.1 Microorganismos 

La cepa axenica numero 15 (posibiemente de Nostoc sp.) y ei consorcio microbiano 
numero 11 (con cianobacterias), fueron proporcionados por ei Laboratorio de Microbioiogia 
Ambientai dei Centro de Investigaciones Cientificas de ia ESPE. 

Los cuitivos fueron mantenidos en medio de cuitivo ifquido BG11 (Rippka, 1988), a un 
fotoperfodo iuz/oscuridad de 16:8, irradiancia de aproximadamente 1000 Lux, a 20±2 ®C de 
temperature y aireacion continue (Koch etal., 2008). 

2.2 Preparacion de agua residual semi-sintetica y sintetica 

Se prepare agua residuai semi-sintetica a partir de agua residuai reai coiectada en ia 
descarga de efiuentes de ia ESPE, sin esteriiizacion previa, supiementandoie fosforo 
inorganico (K2PO3) y carbono organico en forma de peptone y extracto de came (CEE, 1982) 
con ei fin de iiegar a condiciones de un agua residuai urbane media (Metcaif & Eddy Inc. 1995). 

Ei agua residuai sintetica fue preparada de acuerdo a io estabiecido por ia Comunidad 
de Estados Europeos en sus referentes para ios metodos de controi de tensoactivos (CEE, 
1982) asemejando ias concentraciones de nitrogeno amoniacai, fosforo totai y DQO a ios 
parametros de un agua residuai urbane media (Metcaif & Eddy Inc. 1995). 
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2.3 Condiciones del experimento 

Se dispusieron 275 mL de agua residual per unidad experimental. Las condiciones de 
cultivo fueron de 41 61 .25 ± 1 36.37 Lux, a una temperatura de 22 ± 2 -C, con fotoperiodo de 1 6 
horas de luz, 8 de oscuridad (Koch etal., 2008), sin aireacion por 21 dias. Tanto para las 
celulas en suspension e inmovilizadas, se partio de un inoculo inicial de cianobacterias de 
4x1 0^ cel/mL (Chacon et a!., 2004). 

La inmovilizacion celular se llevo a cabo homogenizando los inoculos previamente 
cuantificados, en una solucion de alginato de sodio al 1 ,25% (p/v) mediante goteo sobre una 
solucion de cloruro de calcio al 1% (p/v) con una bureta (Morales, 2006). 

2.4 Diseho experimental 



Se llevo a cabo un DCA de 5 tratamientos con 3 repeticiones para cada uno y pruebas 
de comparacion de mediae de Duncan en el programa SPSS 15.0. 

2.5 Cuantificacion de nitrogeno amoniacal, fosforo total y DQO 

Se utilizaron los metodos colorimetricos de la sal de fenol, el acido 
vanadomolibdofosforico y de reflujo cerrado para la cuantificacion de nitrogeno amoniacal, 
fosforo total y DQO respectivamente (APHA, AWWA & WEF, 1998; APHA, AWWA & WPCF, 
1992). 

2.6 Ensayos de toxicidad en Daphnia pulex 

Se expuso a neonatos de Daphnia pulex con menos de 24 horas de nacidos a cinco 
concentraciones de agua residual semi-sintetica tratada, las cuales fueron de 100%, 50%, 25%, 
12.5% y 6.25%, ademas de un bianco constituido solamente de agua dura como control 
negativo (Castillo, 2004). 



3. RESULTADOS 

3.1 Reduccion de nitrogeno amoniacal, fosforo total y DQO en agua residual 
semi-sintetica 

La concentracion de nitrogeno, descendio de 24.9310.9 mg/L a valores entre 4.5410.9 
mg/L y 0.1610.1 mg/L. En cuanto al fosforo, la concentracion inicial de 11.8710.5 mg/L se 
redujo en todos los tratamientos a valores entre 0.9310.1 mg/L y 0.4010.1 mg/L. Finalmente la 
demanda qufmica de oxfgeno (DQO) se redujo de una concentracion inicial de 411.9313.8 
mg/L a valores entre 1 16.34119.1 mg/L y 60.9614.7 mg/L (Tablal). 

Tabla 1. Concentracion final y porcentaje de reduccion de nitrogeno, fosforo y DQO en agua 
residual semi-sintetica utilizando un cultivo axenico de cianobacterias y un consorcio 
microbiano, inmovilizados y en suspension. 



Tratamiento 


Concentracion final 
mg/L 


Reduccion 

% 


Nitrogeno 


Fosforo 


DQO 


Nitrogeno 


Fosforo 


DQO 


Agua semi-sintetica 


3.4510.5 


0.8410.1 


73.50117.8 


86.1812.3 


92.9111.4 


82.1614.3 


Axenico 


4.5410.9 


0.9310.1 


62.0119.4 


81.7813.6 


92.1710.2 


84.9512.2 


Consorcio 


2.8610.4 


0.9210.1 


60.9614,7 


88.5211.9 


92.2210.9 


85.2011.1 


Ax. Inmovilizado 


0.1610.1 


0.4010.1 


116.34119.5 


99.3710.6 


96.6511.2 


71.7614.6 


Cons. Inmovilizado 


3.5110.0 


0.8410.2 


113.21120.1 


85.9410.3 


92.9612.0 


72.5214.8 
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Para los tres parametros, se hallo alta diferencia significativa entre los tratamientos. 
Realizadas las pruebas de comparaclon de medlas de Duncan, se determino que para la 
reducclon de nitrogeno (99.37%) y fosforo (96.65%) el mejor tratamiento lo constituyo el cultivo 
axenico inmovilizado, mientras que para la reducclon de DQO los mejores tratamientos fueron 
los cultivos suspendidos tanto de la cepa axenica (84,95%) como del consorcio microbiano 
(85.20%), incluyendo el tratamiento del agua semi-sintetica (82.16%). 

3.2 Reducclon de nitrogeno amoniacal, fosforo total y DQO en agua residual 
sintetica 

La concentracion de nitrogeno descendio de 28.72±2.0 mg/L a 22.66±4.5 mg/L 
solamente para el case del tratamiento agua sintetica, mientras que en el resto, el valor final 
fluctuo entre 41.5±5.90 mg/L y 53.23±5.1 mg/L. En cuanto al fosforo, este disminuyo de una 
concentracion inicial de 15.0110.5 mg/L a valores entre 1.0010.6 mg/L y 13.1110.1 mg/L. Por 
otra parte, la demanda qufmica de oxfgeno (DQO), de una concentracion inicial de 324.7613.3 
mg/L paso a valores de concentracion final entre 66.67110.6 mg/L y 29.0114,8 mg/L, 
excluyendo el tratamiento agua sintetica (Tabla 2). 



Tabla 2. Concentracion final y porcentaje de reducclon de nitrogeno, fosforo y DQO en agua 
residual sintetica utilizando un cultivo axenico de cianobacterias y un consorcio 
microbiano, inmovilizados y en suspension. 



Tratamiento 


Concentracion final 
mg/L 


Reduccion 

% 


Nitrogeno 


Fosforo 


DQO 


Nitrogeno 


Fosforo 


DQO 


Agua sintetica 


22.6614.5 


13.1110.1 


324.67135.0 


21.12115.7 


12.6111.0 


0.03110.8 


Axenico 


53.2315.1 


6.3910.1 


66.67110.6 


- 


57.4111.2 


79.4713.2 


Consorcio 


45.1417.2 


5.6310.7 


29.0114,8 


- 


62.4515.1 


91.0711.5 


Ax. Inmovilizado 


43.9611.8 


3.5610.5 


53.4917.1 


- 


76.2913.4 


83.5312.1 


Cons. Inmovilizado 


41.5415.9 


1 .0010.6 


37.4814.8 


- 


93.3614.2 


88.4611.5 



Se encontro diferencia significativa alta entre los tratamientos. Realizadas las pruebas 
de mediae de Duncan, se establecio que el cultivo del consorcio inmovilizado fue el que mejor 
eficiencia de reducclon de fosforo total presento (93.36%), mientras que su par suspendido, 
mostro la mejor eficiencia de reducclon de DQO (91.07%). Para el case del nitrogeno 
amoniacal, ninguno de los tratamientos, a excepcion del de agua sintetica, demostro reducclon 
alguna con respecto a la concentracion inicial del contaminante. 



3.3 Pruebas de toxicidad en Daphnia pulex 



Transcurridas las 48 horas de exposicion de los organismos de prueba al agua tratada 
correspondiente al tratamiento del cultivo axenico inmovilizado, no se hallo mortalidad de los 
individuos a ninguna de las concentraciones utilizadas 



4. DISCUSION 

Cuando la taza nitrogeno-fosforo esta en un range entre 2 y 5, como en el presente 
trabajo, es de esperarse una predominancia en el crecimiento de las cianobacterias con 
respecto a otros grupos taxonomicos existentes en el medio (Levich, 1996). Por tanto, serfa 
atribuible que la reducclon de estos dos elementos estarfa relacionada estrechamente con el 
crecimiento predominante de cianobacterias. Sin embargo, se debe considerar el rol que 
cumplieron los microorganismos nativos propios del agua residual en los procesos de 
degradacion, como se evidencio en el ensayo con el agua residual sintetica que no presento 
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reduccion de nitrogeno, comparado con el llevado a cabo en el agua residual seml-sintetica, 
fenomeno citado por Moreno etal. (2010, en prensa). La mayor reduccion de este elemento 
fue deblda a la presencia de otros microorganismos, ya que el amonlo representa la fuente de 
nitrogeno mas accesible tanto para autotroficos como para heterotroficos (Singleton, 2004; 
Muro-Pastor et al., 2005; Flores et al., 2005). 

Los cultivos Inmovlllzados, presentaron mayor reduccion de nitrogeno y fosforo en 
relaclon a la observada por sus pares en suspension, de manera similar a lo citado por Mallick 
(2002) en un estudio realizado con la microalga verde Chlorella vulgaris y la clanobacterla 
Anabaena doliolum. 

En el tratamlento con agua sintetica, sin Inoculo microbacterlano alguno, la reduccion 
de nitrogeno se deblo a la volatlllzaclon del nitrogeno por procesos fisico-quimicos de equillbrio, 
al presentarse en el medio lones hidroxllo por un pH llgeramente basico (Loomis & Connor, 
2002 ) 

Las clanobacterlas, a pesar de ser microorganismos fotoautotroficos, son capaces de 
utilizer carbono organico en sus procesos metabollcos (Yu, JIa & Dal, 2009; MIshra, Kumar & 
Jha, 1985) lo que explIca la reduccion considerable de DQO presentada por el cultivo axenico 
suspendido en el agua residual sintetica. 

La ausencia de toxicldad del agua tratada sobre Daphnia pu/expuede expllcarse a, que 
como se ha descrito, (Leflalve & Ten-Hage. 2007) las toxinas producldas por las clanobacterlas 
de la famllla Nostocales son llberadas al medio cuando se produce llsis celular, como resultado 
de crecimlento masivo, lo cual no se presento en la presente Investigaclon, o a que esta cepa 
de clanobacterlas no presentan toxicldad. 

5. CONCLUSIONES 

El cultivo axenico de clanobacterlas Inmovlllzado presento la mayor eficlencia de 
remoclon de nitrogeno amonlacal y fosforo total, mientras que los tratamlentos en suspension, 
asf como el agua residual seml-sintetica con la microblota original del agua residual, 
presentaron una capacldad similar de reduccion de DQO. 

El ensayo en agua residual sintetica demuestra que el consorcio microblano y la cepa 
de Nostoc sp. son capaces de reduclr las concentraclones de fosforo y DQO 
Independlentemente de la microblota propla del agua residual. 

El agua residual tratada medlante un cultivo axenico de Nostoc sp. Inmovlllzado no 
presento toxicldad en Daphnia pulex. 
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RESUMEN 



Se describe con cierto detalle el puente Carrizal, ubicado en la Provincia de Manabi, 
que esta compuesto por dos estribos abiertos con contrafuerte, uno de ellos se encuentra 
sobre pilotes. Adicionalmente existen dos apoyos interiores, cada uno de ellos conformado por 
9 pilotes circulares y cuatro niveles de vigas de arriostramiento dos de ellas son de acero y dos 
de hormigon armado. 

Posteriormente se describen varies programas que se ban implementado para el 
software CEINCI-LAB. Estos programas son para el calculo de las matrices de rigidez y de 
masas del estribo izquierdo, de uno de los apoyos, del estribo derecho y de una de las vigas de 
acero de 75 m., de longitud. En el anexo 1 se indica el programa final con el cual se realizo el 
analisis sfsmico longitudinal empleando los programas de CEINCI-LAB. 

Se comparo la respuesta en el tiempo, empleando dos modelos de analisis, en el 
primero se considero el empuje active del suelo sobre los estribos y en el segundo modelo no 
se considero este empuje active. Para el primer modelo el vector de condiciones iniciales no es 
nulo y esta compuesto por los desplazamientos laterales que genera el suelo en el estribo. 
Para el segundo modelo el vector de condiciones iniciales es nulo. Se utilizo un sismo artificial 
de 20 s., de duracion compatible con el espectro del Codigo Ecuatoriano de la Construccion 
para suelo S3. 



ABSTRACT 

It describes in some detail the bridge Carrizal, located in the province of Manabf, which 
consists of two open abutments with buttresses, one of them is on piles. Additionally two interior 
supports exist, each one of them conformed by 9 circular piles and four levels of bracing beams 
two of them is of steel and two of armed concrete. 

Subsequently described several programs that have been implemented to CEINCI-LAB 
software. These programs are for the calculation of stiffness matrices and mass of the left 
abutment of one of the supports, the right abutment and one of the beams steel of 75 m. long. 
Annex 1 show the final program which was conducted using longitudinal seismic analysis 
programs CEINCI-LAB. 

We compared the response in time, using two analysis models, the first was considered 
active soil thrust the abutment and the second model did not consider this passive resistance. 
For the first model, the vector of initial conditions is not null and is composed of the lateral 
displacements generated by the soil in the abutment. For the second model, the vector of initial 
conditions is null. We used artificial earthquake 20 seconds, time compatible with the spectrum 
of the Ecuadorian Code of Construction for soil S3. 
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11.DESCRIPCION DE PUENTE 



En la Provincia de Manabi se encuentra el Puente Carrizal, construldo per el Cuerpo de 
Ingenleros del Ejercito y forma parte de la carretera QuIroga-PIchIncha. En la parte superior de 
la figura 1 se apreclan las poblaclones de Quiroga y PIchIncha y la via se encuentra dentro del 
circulo, al lado del campamento se halla ubicado el puente. 

El puente tiene una longitud de 249.20 m., esta compuesto por dos estrlbos y dos pllas 
Intermedlas, como se observe en la parte Inferior de la figura 1. En los espaclos llbres de cada 
tramo se tienen cuatro vigas metallcas de 75 m., de longitud, unidas por diafragmas y tensores 
Inferlores. Las pllas Intermedlas estan sobre nueve pllotes prebarrenados de 2.0 m., de 
diametro que cubren tramos de 12 m., del tablero del puente; en cada tramo esta previsto la 
construcclon de juntas de dllataclon de 5 cm., con lo que se tiene los 249.20 m. Morales (2009). 
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Figura 1 Panoramica del puente sobre el rfo Carrizal. (Cortesfa de E. Morales 2009). 



En la figura 2 se apreclan los dos estrlbos del puente; el que se encuentra en la parte 
superior no tiene pllotes de cimentacion; se aprecia que en la montafia se ha realizado un corte 
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para el eje de la via (fotografia superior Izquierda), a este estrlbo denomlnaron Estribo Quiroga, 
ya que se encuentra mas cerca de la poblaclon de Quiroga. La parte central del estrlbo es recta 
y tiene dos alas en sus lados. En la fotografia superior derecha se presenta otra vista del 
estrlbo, en la cual aparecen las trabas sfsmicas, en la parte superior, que Impiden que las vigas 
de acero de 75 m., de longitud se salgan ante movimlentos en sentido transversal al eje del 
puente, ya que las vigas de acero se encuentran simplemente apoyadas sobre las vigas del 
estrlbo, tema que se detallara mas adelante. 

En la parte Inferior de la figura 2, aparece el otro estrlbo que si tiene pllotes en su 
cimentaclon. En la fotograffa Inferior Izquierda se observan los cabezales de los pllotes que 
fueron fundidos conjuntamente con la viga cabezal, en la fotograffa aparecen 14 pllotes no se 
apreclan 2 mas; en total son 16 pllotes. En la fotograffa Inferior Izquierda se aprecia las vigas 
metallcas que se ban colocado sobre el estrlbo; a este lo llaman Estribo Pichincha. En esta 
ultima figura se observe la construcclon de los dos apoyos Interlores del puente. 





Figura 2 Estribos del puente sobre el rfo Carrizal. (Cortesfa de Morales E., 2010). 



La geometrfa de los dos estribos es practicamente la misma por eso se describe en 
forma general los elementos estructurales que contlene cada uno de ellos. Como se manifesto, 
el estrlbo tiene una parte central recta y dos alas, como se observe a la Izquierda de la figura 3, 
que corresponde a una vista posterior del estrlbo; se observe que hay dos contrafuertes en 
cada una de las alas y cuatro contrafuertes en la parte central; en total se tienen echo 
contrafuertes de 0.40 m., de espesor con una base Inferior de 2.85 m., y una base superior de 
0.85 m.; la altura de estos contrafuertes es de 5.93 m. Para el anallsis sfsmico que se presenta 
en este artfculo, los contrafuertes fueron modelados como rectangulares. VInueza (1992) 
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Figura 3 Vista posterior y iateral de ios estribos. 



En la parte frontal de Ios contrafuertes se encuentra una pantalla de 0.35 m., de 
espesor y de la misma altura de Ios contrafuertes; sobre estos dos elementos estructurales 
(contrafuerte y pantalla) hay una viga de 1 .20 m., de base y 0.60 m., de alto. Es sobre esta viga 
que se van a asentar Ios apoyos de neopreno y las cuatro vigas de acero tlpo I , que como se 
Indico, tienen una longitud de 75.0 m., cada una. 

Para termlnar con la descripcion rapida del estrlbo se debe IndIcar que en la parte 
superior hay otra pantalla de 0.30 m., de espesor y 3.30 m. de alto en la parte recta ya que en 
las alas, la altura es variable como se ve en la figura 3. Por ultimo pero en primer lugar se debe 
IndIcar que la cimentaclon de la pantalla tiene 1.20 m., de alto en el estrlbo con pllotes y 1.00 
m., de alto en el estrlbo sin pllotes; en la figura 3 se aprecia la forma de la cimentaclon que 
dicho sea de paso se encuentra sobre un replantlllo de 0.10 m. 

Una vez que se ha descrito la geometrfa de Ios estribos, se pasa ahora a IndIcar la 
forma de Ios apoyos centrales del puente. Se habfa Indicado que cada apoyo tiene nueve 
pllotes, Ios mismos que se presentan en la parte superior de la figura 4, a la derecha de esta 
figura se observe las camisas de acero de Ios pllotes que tienen un espesor de 1 .5 cm., fueron 
colocadas las camisas con el objeto de protegerlas del medio amblente pero estas camisas 
proporclonan rigidez y resistencia adiclonal a la del pllote de hormigon armado, de 2 m., de 
diametro. 

En la fotograffa Inferior de la figura 4 se observe Ios pllotes y las vigas de 
arriostramlento superior e Inferior de acero, que se hallan en las cotas 56 m., y 62 m., las cotas 
menclonadas son referldas a la parte Inferior de la viga; estas vigas tienen un peralte de 1 .60 
m., y 0.90 m., respectivamente, son vigas tlpo I (la cota de desplante de Ios pllotes es 16.82 
m. Morales 2009). 

De Ios nueve pllotes que existen en cada apoyo, sels son exterlores y tres Interlores. 
Ahora blen, Ios pllotes exterlores son de 2 m., de diametro en una altura de 54 m., y luego se 
reducen a 1 m., de diametro como se observe en la figura 5. Los tres pllotes Interlores son de 2 
m., de diametro en toda su longitud que esta alrededor de Ios 58.50 m. 

Dos aspectos son Importantes de destacar con relaclon a la figura 5. El primero que a 
mas de las vigas de arriostramlento de acero, existe en Ios apoyos dos vigas mas de hormigon 
armado, la una se halla en la cota 70.279; es una viga de 0.80 m., de base por 1.25 m., de 
altura; la otra viga tiene una base de 0.50 m., y una altura de 0.80 m. El segundo aspecto es 
que la reducclon del diametro de Ios pllotes exterlores de 2.0 m., a 1 m., es para tener especio 
para apoyar las vigas longitudinales de acero de 75 m., por este motivo en la figura 5 se 
observe que la viga de hormigon armado de 1 .25 m., de peralte tiene unos voladizos para que 
el apoyo tenga una longitud de 1 .0 m. 





Figura 4 Pilotes y vigas de acero de apoyos centrales. 
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En la figura 5 se ha colocado el apoyo de neoprene y la viga de acero de 75 m., de 
longitud unicamente en el lado derecho pero tamblen existe lo mismo en el apoyo Izquierdo, no 
se lo dibujo para no compllcar la expllcaclon. Per este ultimo motive no se Indica nada de las 
vigas transversales, que se apreclan en la figura 5. 



12.PROGRAMAS ESPECIFICOS DE CEINCI-LAB 



CEINCI-LAB es un software educative y practice desarrollado por Aguiar (2010) que 
contlene una serle de programas para el: anallsis estatico de estructuras, anallsis sfsmico de 
puentes, estructuras con alsladores de base, entre otros temas. 

Con el desarrollo de este artfculo se Incorporan por una parte nuevos programas a 
CEINCI-LAB, los mismos que se IndIcan en este apartado y por otra parte se realize el anallsis 
sfsmico complete en sentido longitudinal en los sigulentes apartados, utlllzando todos los 
programas existentes en la carpeta denomlnada Carrizal de CEINCI-LAB 



12.1.1 Estribo Izquierdo 



En la figura 6, se muestra el modelo numerico de anallsis del estribo Quiroga que se halla 
sobre 16 pllotes, cuatro en el ala Izquierda, echo en la parte central y cuatro en el ala derecha. 
Se realize un anallsis sfsmico en 2 dlmenslones, por esa razon es que se apreclan dos pllotes. 




Air 









11 



—13 @45 cm — 
5,85 



Figura 6 Modelo numerico de estribo izquierdo. 

Lo que interesa ilustrar es que los pllotes han side modelados con elementos finitos 
lineales que tienen tres grades de libertad por nudo, un desplazamiento horizontal, un 
desplazamiento vertical y una rotacion. Se consideraron 10 elementos finitos en cada pilote; se 
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obtuvo la matriz de rigidez de los pilotes de 60 per 60 y luego se condense a 20 por 20. Aguiar 
(2004). Los grades de libertad que se indican a la izquierda de la figura 6, referentes a los 
pilotes, son las coordenadas principales. La matriz de masa de los pilotes es de 20 por 20 y se 
halla en forma similar al de la matriz de rigidez pero utilizando la matriz de masa del elemento 
finite de masa uniforme distribuida. 

En el estribo, propiamente dicho, se trabajo con elementos finitos rectangulares con 
dos grades de libertad por nudo. La base de los estribos se halla apoyada sobre el suelo y otra 
sobre los pilotes. En este case las masas se concentre en puntos discretes sehalados con un 
cfrculo. 



Para el estribo se indica a la izquierda de la figura 6 todos los grades de libertad 
contados a partir de la ultima coordenada principal de los pilotes. Se hallo la matriz de rigidez 
de toda la estructura de 168 por 168 y se condense a las 44 grades de libertad indicadas a la 
derecha de la figura 6. 

El programa que obtiene la matriz de rigidez y de masas (de 44 por 44) en las 
coordenadas principales indicadas a la derecha de la figura 6 se denomina 

rigidez_estribo_izquierdo. 



12.1.2 Apoyos Centrales 



En dos dimensiones se aprecian tres pilotes y cuatro vigas; las vigas superiores son de 
hormigon armado y las inferiores son de acero. De paso se indica la forma como se numeran 
los grades de libertad en CEINCI-LAB, para el case de los pilotes. Todo esto lo hace el 
programa, solamente se indica para que el lector comprenda como se trabajan los pilotes. A la 
izquierda de la figura 7 se muestran las coordenadas totales de la estructura y a la derecha las 
coordenadas principales. 

Primero se numeran todos los desplazamientos horizontales, empezando desde el 
pilote izquierdo, de abajo hacia arriba; despues todos los grados verticales y finalmente las 
rotaciones. Las vigas del apoyo central, tambien han sido analizadas con elementos finitos 
lineales con masa uniforme distribuida. 

Para usar los programas de la librerfa de CEINCI-LAB, primero se numeran los 
elementos de los pilotes; luego se numeran las vigas empezando desde la viga superior como 
se observe en la figura 8; los numeros que estan fuera del cfrculo corresponden a la 
numeracion de los nudos y los que se encuentran dentro del cfrculo a la numeracion de los 
elementos. El programa que determine la matriz de rigidez y la matriz de masas del apoyo 
central, en CEINCI-LAB se denomina rigidez_apoyo_central 

Si el usuario de CEINCI-LAB quiere trabajar con mas elementos finitos en los pilotes. 
Tiene que hacer un nuevo programa pero basado en rigidez_apoyo_central , para lo cual se 
recomienda cambiar de nombre el programa e indicar el numero de elementos finitos {ne). 
Pero tendra que dibujar la estructura para saber cual es el nudo inicial y el nudo final de las 
vigas, informacion que debe modificar al programa inicial en el cual esta trabajando. No hay 
mas cambios. 

Si tiene dos ejes de pilotes o mas de tres ejes de pilotes, en el nuevo programa que 
esta realizando, utilizando rigidez_apoyo_central, debe indicar el numero de ejes de pilotes 
que tenga (np ) ; Encontrar una ecuacion para encontrar el numero total de elementos de la 
estructura y cambiarla por que la que esta actualmente es para el case de 44 elementos. Hay 
mas cambios que se deben realizar para ello se recomienda que se vea todo el programa y se 
vayan realizando los cambios respectivos. Este al igual que todos los programas de CEINCI- 
LAB, estan muy bien comentados. 
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Figura 7 Coordenadas Totales y Coordenadas principales de un apoyo central. 
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Figura 8 Numeracion de nudes y elementos de apoyo central. 
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Los elementos del 1 al 11 del apoyo central tienen un diametro de 2.0 m., el elemento 
12 tiene un diametro de 1.0 m., por que hay una reducclon de la secclon para que se apoyen 
las vigas de 75 m. Todos los elementos de la plla central son de 2.0 m., en la tercera flla los 
elementos que van desde el 25 al 35 son de 2.0 m., de diametro y el elemento 36 de 1 .0 m., de 
diametro. En el programa rigidez_apoyo_central , se programo esto en pocas lineas, 
basicamente lo que se hizo fue un bucle desde el elemento 1 hasta el elemento np*ne y 
dentro de este bucle hay una pregunta iEs el elemento 12 o 36?, en case afirmativo el 
diametro es 1 .0 m., caso contrarlo 2.0 m. 





Figura 9 Modelo numerico de estribo derecho. 

Para las vigas de hormigon se IndIcan la base y la altura de cada una de ellas. Para las 
vigas de acero se Indica el area y el memento de Inercia de la secclon transversal. Hay varlas 
formas de Ingresar estos dates, una de ellas es la que consta en rigidez_apoyo_central. 

12.1.3 Estribo derecho 

El estribo derecho es similar al estribo Izquierdo pero dado la vuelta y la cimentaclon no 
tiene pllotes. Todo lo Indicado en 2.1.1 es apllcable para este estribo, por lo que se van a dar 
algunos comentarlos adiclonales de como se realize el programa rigidez_estribo_derecho. 

La primera flla de elementos finitos corresponde a la cimentaclon que tiene una altura 
de 1.0 m., la profundidad de este elemento que no se ve en la figura 9 es de 12.0 m. Para la 
cimentaclon se pudo haber trabajado con 2 elementos fInItos en sentido Y, con lo que la altura 
de cada elemento en sentido vertical serfa de 0.50 m. 

Para los elementos finitos que se hallan sobre la cimentaclon se deben agrupar los 
elementos finitos que corresponden a la pantalla y los elementos finitos que son de los 
contrafuertes; debldo a que la profundidad es diferente. Para los elementos finitos de la 
pantalla esta profundidad es de 11.0 m.; en camblo para los contrafuertes la profundidad es 
1.60 m. Notese que en la profundidad se han considerado 4 contrafuertes que son los que se 
encuentran en la parte central del estribo. (HIpotesIs de calculo). 



Revista CIENCIA Volumen 13 N.- 1 de 2010 



71 



La profundidad de los contrafuertes se obtuvo sumando su espesor (Criterio adoptado 
para encontrar el espesor del conjunto). Se puede apllcar otros criterlos mas refinados para 
encontrar el espesor equivalente uno de ellos es en base al memento de Inercla. En este case 
primero se debe determiner el centre de gravedad de todos los contrafuertes a una 
determlnada altura ya que son de secclon variable; segundo se halla el memento de Inercla de 
todos los contrafuertes con respecto al centre de gravedad y finalmente se halla el espesor 
equivalente de tal manera que se tenga el mismo memento de Inercla. Aguiar y Espinoza 
( 2010 ). 



Luego se dan los dates de la viga superior sobre la cual se aslentan las vigas de acero 
de 75 m., y que soportan a la pantalla superior. La altura de esta vIga es de 0.60 m. SI se 
desea se puede considerar 2 elementos finitos en sentido vertical en lugar de uno, como se 
calculo; en este case la altura de cada elemento es 0.30. Finalmente, se dan los datos de los 
elementos fInItos de la pantalla superior. 

12.1.4 Vigas de acero 

En cada viga de acero se considero 150 divislones en sentido X (Luz de 75 m.) y 3 
divislones en sentido Y (Altura de la viga de 3.0 m.). De tal manera que en cada viga se ha 
trabajado con 450 elementos finitos de 0.50 m., por 1.00 m. En la figura 10 se muestra la 
numeraclon de los grados de llbertad, las masas se concentraron en los nudos exterlores por 
este motive se numeran primero los grados de llbertad de estos nudos, ya que estas son las 
coordenadas principales. Notese que en sentido horizontal se tienen 8 coordenadas principales 
y en sentido vertical se tienen 6 coordenadas principales mas. Luego se numeran todos los 
demas grados de llbertad. En la parte Inferior de la figura 10 se Indica un detalle del apoyo de 
Neopreno, (Toapanta 2009), el mismo que ha sido modelado como un resorte de rigidez lineal. 

SI se desea realizar el anallsis sfsmico para la componente horizontal y vertical de 
movimlento del suelo, una vez obtenida la matriz de rIgIdez de la estructura por ensamblaje 
directo se condense a una matriz de 14 por 14. Ahora, cuando solo se realize el anallsis 
sfsmico para la componente horizontal la condensaclon se hace a una de 8 por 8. 

En la figura 1 1 se IndIcan las secclones transversales del perfll tlpo 7 a lo largo de la 
viga, en los extremos tiene una menor secclon y en el tramo central una mayor secclon. En la 
parte Inferior de la grafica se muestra la abscise X., hasta el centre de luz ya que es simetrica 
la viga. En la table 1 se IndIcan los espesores equivalentes como si se tratara de una viga 
rectangular que tiene el mismo memento de Inercla de la viga I . 





Figura 10 Numeraclon de los grados de llbertad de una viga. 
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Figura 1 1 Secciones transversales de las vigas de acero 




DeX= 13.31 
DeX= 15.81 



De X= 15,81 
De X= 37,50 



Tabla 1 Espesores equivalentes. 



JC 


0 a 0.5 m. 


0.5 a 6.9 m. 


6.9 a 10.81 
m. 


10.81 a 13.31 
m. 


13.31 a 15.81 
m. 


15.81 a 37.50 
m. 


e (cm.) 


3.71 


3.09 


3.57 


4.48 


4.83 


4.99 



El programa que determina la matriz de rigidez y de masas en coordenadas principales 
se denomina: rigidez_viga_superior. 



13.MODELO NUMERICO DE ANALISIS 



En la figura 12 se muestra una vista en ei sentido de anaiisis sismico, que se considera 
en este articuio. Notese que ei estribo izquierdo tiene piiotes y ei estribo derecho no ios tiene. 
De iguai manera se aprecia que ios piiotes de ios apoyos centraies tienen iongitudes mucho 
mas aitas que ios piiotes dei estribo izquierdo. 




En ia parte superior de ia figura 13 se muestran ias coordenadas principaies (grades 
de iibertad horizontai) dei puente Carrizai. Para que se pueda ver de mejor forma, en ia parte 
superior se indican ios grades de iibertad dei estribo izquierdo con sus respectivos piiotes, de 
ios apoyos centraies y dei estribo derecho. 

En ia parte inferior de ia figura 13, se muestran ios grades de iibertad horizontaies de 
ias vigas de acero; ia primera corresponde a ia viga que se encuentra entre ei estribo izquierdo 
y ei apoyo centrai; ia segunda a ia viga que se haiia entre ios dos apoyos y ia tercera a ia viga 
que se haiia entre ei segundo apoyo y ei estribo derecho. Es importante saber identificar ios 
grades de iibertad, considerados para ei reporte de resuitados que se daran posteriormente. 

En totai se tiene 164 coordenadas principaies. Luego ias matrices de rigidez, masa y 
amortiguamiento con ias cuaies se reaiiza ei anaiisis sfsmico son de ese orden. Se debe 
destacar que ei sueio fue modeiado mediante resortes de comportamiento no iineai. Las 
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propiedades de la curva carga deformacion del suelo estan Indicadas en Aguiar et al. (2010) y 
fueron Inferldas del estudio de suelos realizado per Beltran (2009), Se considero que el suelo 
tiene un comportamlento elasto perfectamente plastico. 





APOYO CENTRAL 1 APOYO CENTRAL 2 





Figura 13 Coordenadas Principales de modelo de analisis sismico longitudinal. 
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14.EMPUJE ACTIVO DEL SUELO 



Cuando la estructura esta en reposo, existe un empuje del suelo sobre los dos estrlbos 
del puente; empuje que genera desplazamlentos en estos elementos estructurales. En el 
articulo se encontro la respuesta sismica de dos formas, considerando la preslon que el suelo 
transmite a los estrlbos y sin considerar esta preslon. 

Para el primer modelo se hallo la fuerza horizontal F en cada uno de los grados de 
llbertad del estrlbo Izquierdo, como se Indica en la parte Izquierda de la figura 14. Luego se 
obtuvo el vector de cargas generallzadas Q y se calculo los desplazamlentos horizontales q 
en dichos grados de llbertad. Aguiar (2004). 



F = 



"T" 



bs 



(1 ) 



Donde y es el peso especffico del suelo; h es la altura desde la superficle al punto 
donde se desea calcular la fuerza F ■, es el ancho del volumen de suelo considerado en el 
anallsis, se trabajo con = 10 m. 




Figura 14 Fuerzas estaticas actuando en los estrlbos, debido al suelo. 



Cuando las cargas actuan en los nudos, la obtenclon del vector de cargas es directo. 
Ahora blen, en la figura 6 se muestran las coordenadas principales del estrlbo Izquierdo. Por lo 
tanto, se puede trabajar con un vector de cargas Q de 44 elementos y hallar los 
desplazamlentos horizontales q , resolvlendo el sigulente sistema de ecuaclones llneales. 




Q = KE I q 



( 2 ) 



Slendo KEI la matriz de rigidez del estrlbo Izquierdo condensada a las coordenadas 
laterales de 44 por 44. Q el vector de cargas compuesta por valores en los grados de llbertad 
donde empuja el suelo y los demas nulos; q es el vector de desplazamlentos de 44 por 1 . 
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En forma similar se trabaja en el estribo derecho. En este caso las fuerzas F son 
negativas debido a que las fuerzas del suelo estan actuando en sentido contrario a la direccion 
positiva de los grades de libertad. Se hallan los desplazamientos en el estribo derecho y con 
todos los desplazamientos se tiene el vector de condiciones iniciales q para t = 0 . 

Por lo tanto, en el modelo 1 , se obtuvo la respuesta en el tiempo con condiciones 
iniciales diferentes de cero; es para el caso en que se considera empuje de suelo y en el 
modelo 2, se encontro la respuesta con vector ^ = 0 para t = 0 , es para el caso en que no se 
considera el empuje de suelo. 



15.ELEMENTOS FINITOS Y SISMO DE ANALISIS 



Los pilotes fueron modelados con elementos finitos lineales, con masa uniforme 
distribuida. Aguiar (1981). Las funciones de forma con las cuales se halla la matriz de rigidez y 
de masas, son las siguientes: 












= (/>,(x) = l-3 — +2— <!>,{x) = X 






X‘ 






3 - 2^1 

L) 



1-^ 

Lj 






X 



( 3 ) 



1 

Lj 



SiendoL, la longitud del elemento y X la abscisa de un punto del elemento. Los elementos 
de la matriz de rigidez A: , y de masas m de un elemento, se obtienen a partir de las funciones 
de forma, con las expresiones siguientes: 

L L 

k(i, j) = ^E1 (j). ^jdx m(i, y) = J m (j). dx ( 4 ) 

0 0 



Donde El es la rigidez a flexion del elemento; m es la masa por unidad de longitud. 
Se obtuvo la matriz de rigidez y de masas por ensamblaje directo y luego se condense estas 
matrices a los grades de libertad horizontales. Aguiar (2004, 2007). 

Para los estribos propiamente dicho se trabajo con elementos finitos rectangulares; 
con dos grades de libertad por nudo, desplazamiento horizontal y vertical; con el proposito de 
contemplar la flexion en el centre de gravedad del elemento se considero cuatro grades de 
libertad mas. Las funciones de forma utilizadas en coordenadas naturales se presenta en la 
tabla 2. Ohate (1992). 

Se trabajo con dos puntos de la cuadratura de gauss, en cada sentido, para hallar la 
matriz de rigidez del elemento finite, con la siguiente ecuacion: 



k = \B‘EBdV (5) 

V 

Donde k es la matriz de rigidez de un elemento finite; V es el volumen; B es la 
matriz que relaciona los esfuerzos con las deformaciones es una matriz de 3X12; E es la 
matriz de elasticidad para el caso de tension plana. En funcion del espesor e del elemento, la 
ecuacion ( 5 ) se transforma en: 
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k = E Bdia[j)edA (6) 

- 1-1 

Las variables no definidas son dA que es el diferencial de area; 7 la matriz 
Jacobiana; det(7) es el determinante de la matriz Jacobiana, que reporta el area de influencia 
en el punto de la cuadratura de gauss. 



Tabla 2 - Funciones de forma para elementos finitos rectangulares 



Coordenada natural 


Coordenada 

Centroidal 


Deformada 

elemental 


Funcion de Forma 


1 




ul 


=0.25(-5 + l)(-r + l) 


vl 


=0.25(-5 + l)(-r + l) 


2 


u2 


= 0 . 25(5 + l)(-r + l) 


v2 


= 0 . 25(5 + l)(-r + l) 


3 


u3 


^3 = 0 . 25(5 + l)(r + l) 


v3 


(^3 = 0 . 25(5 + l)(r + l) 


4 


u4 


(z)4=0.25(-5 + l)(r + l) 


v4 


(^4=0.25(-5 + 1)(? + 1) 




5 


u5 


(^ 5 = 1 - 5 ^ 


v5 


(Z>5 = 1 - 5 ^ 


6 


u6 




v6 





En la figura 15 se indica el acelerograma artificial utilizado, tiene una duracion de 20 s., 
una aceleracion maxima del suelo de 0.4 g ., siendo g la aceleracion de la gravedad. A la 

derecha de la figura 15 se muestra el espectro de respuesta elastico para un 5% de 
amortiguamiento. El acelerograma artificial generado es compatible con el espectro del Codigo 
Ecuatoriano de la Construccion, CEC-2000 para un perfil de suelo S3. 




Figura 15 Acelerograma utilizado y espectro elastico. 



16.RESULTADOS EN ESTRIBO IZQUIERDO 

En la parte superior de la figura 16 se muestra la historia de desplazamientos 
horizontales en el estribo izquierdo, hallados con el modelo 1 y en la parte inferior de dicha 
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figura la respuesta en el tiempo encontrada con el modelo 2. Como era de esperarse son 
mayores los desplazamlentos encontrados con el modelo 1 (empuje activo del suelo). En 
termlnos generales se puede IndIcar que el estrlbo es bastante rigido y los desplazamlentos 
esperados son bastante bajos. 






Figura 16 Respuesta en Estrlbo Izquierdo con dos modelos. 



Los grados de llbertad selecclonados son el 10 (Punto superior de los pllotes), 37 (Bajo 
el apoyo de Neopreno), 141 (Sobre el apoyo de Neopreno), 147 (En la parte superior de la viga 
de acero, del apoyo de Neopreno). Es Interesante notar la diferencia de desplazamlentos bajo 
el apoyo de Neopreno y sobre el apoyo de Neopreno; los desplazamlentos sobre el apoyo son 
mayores y es debldo a que el Neopreno trabaja como un alslador elastomerico. 

En la figura 17 se muestran los resultados encontrados con el modelo 2 (no considera 
empuje del suelo) en los mismos grados de llbertad. Notese que las fuerzas horizontales son 
bastantes bajas, debldo a que los desplazamlentos laterales son pequefios. 
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Figura 17 Fuerzas horizontales encontradas con modelo 2. 



17.RESULTADOS EN APOYOS CENTRALES 



En la figura 1 8 se muestran la respuesta en el tiempo de desplazamientos hallados con 
el modelo 1, en el apoyo central. Se presenta en los grades de libertad 52 (viga de 
arriostramiento inferior de acero), 53 (viga de arriostramiento superior de acero); 55 (viga 
inferior de hormigon armado) y 56 (viga superior de hormigon armado). 





Figura 18 Respuesta de desplazamientos en apoyo central. 



Los desplazamientos hallados con el modelo 2 son similares a los encontrados con el 
modelo 1 , razon por la que se omite su presentacion. El desplazamiento maximo en el grade de 
libertad 52 esta alrededor de los 25 cm.; el grade de libertad 53 esta en el orden de los 30 cm., 
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el 55 en 35 cm., y el 56 en 36 cm. Indudablemente que la presencia de las vigas ha reducldo 
los desplazamlentos relatives especlalmente entre los grades de llbertad 55 y 56. 

Ceme era de esperarse les desplazamlentes laterales en el apeye central sen muche 
mas altes que les desplazamlentes laterales de les estrlbes. Es Impertante destacar que en les 
des medeles ne se ha Inclulde el amertiguamlente que preperclena el Neeprene a la estructura, 
unicamente se le medele ceme un reserte cuya rigidez es funclon del area de la gema que 
trabaja al certe, del module de corte de la gema y de la sumatorla de la altura de las gomas. Al 
considerar el amertiguamlente en el medele de seguro que van a disminuir los desplazamlentos 
laterales. 




0 5 10 15 20 0 5 10 15 20 




Figura 19 Fuerzas horizontales en apoyo central. 



En la figura 19 se muestra la respuesta de las fuerzas horizontales, a nivel de las vigas 
de acero son bastante altas y a nivel de las vigas de hormigon disminuyen, especlalmente en la 
viga superior. 



18.RESULTADOS EN ESTRIBO DERECHO 



Los resultados encontrados en el estrlbo derecho, son simllares a los hallados en el 
estrlbo Izquierdo, en el sentido de que con el modelo 2 los desplazamlentos son menores en 
relaclon al modelo 1 y que estos desplazamlentos son mucho menores a los desplazamlentos 
encontrados en los apoyos. 

En la figura 20 se muestran los resultados de desplazamlentos, los Indicados en la 
parte superior corresponden al modelo 1 y los Indicados en la parte Inferior al modelo 2. Se 
presenta para el grade de llbertad 118 (Inicio de la pantalla Inferior del estrlbo); 132 (bajo el 
apoyo de Neopreno); 158 (sobre el apoyo de neopreno) y 164 (sobre la vIga de acero de 75 m., 
de longitud). 

Las fuerzas horizontales en el estrlbo derecho son menores a las encontradas en el 
estrlbo Izquierdo y con un comportamlento muy similar. 
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Figura 20 Desplazamientos hallados con los dos modelos en estribo derecho. 



19.CONCLUSIONES 



Se inicio el articulo detallando la forma de los estrlbos del puente sobre el rio Carrizal, 
tanto del que se encuentra sobre pllotes como el que no lo esta. Posterlormente se paso a 
describir la geometria de los apoyos del puente compuesta por nueve pllotes circulares y cuatro 
fllas de vigas, finalmente se expllco la geometria de las vigas de acero de 75 m., de longitud 
que se sustentan en los estrlbos y los apoyos del puente de 292.20 m., de longitud. Era 
necesarla esta descripcion para entender los modelos numericos de anallsis, tanto a nivel de 
subestructuras como de estructura completa. 

Para cada una de las subestructuras se ha presentado el modelo numerico con el cual 
se obtuvo un programa de computaclon que ha sido Incorporado el programa CEINCI-LAB que 
permite encontrar las matrices de rigidez y de masas en coordenadas laterales. Al final en el 
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anexo 1, de este articulo se muestra el programa complete para el anallsis sismico del puente 
Carrizal, utlllzando la llbreria de programas de CEINCI-LAB. Se destaca que todos los 
programas para el anallsis sismIco de un determlnado puente se hallan en una carpeta con el 
nombre del puente. 

FInalmente, se comparo la respuesta en termlnos de desplazamlentos y de fuerzas de 
la respuesta que se obtuvo al realizar el anallsis sismico conslderando el empuje lateral del 
suelo sobre los estrlbos (modelo 1) y al no considerar el empuje lateral del suelo (modelo 2). 
Evidentemente que los desplazamlentos obtenidos con el modelo 2 fueron llgeramente 
mayores a los encontrados con el modelo 1 , en los estrlbos. Por lo tanto, el empuje active del 
suelo sobre los estrlbos debe ser modelado en todo anallsis sismico de puentes. 
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ANEXO 1 



clear all; close all; echo off 

% 



% Programa para el Anallsis Sismico Longitudinal de Puente Carrizal 

% 

% Viga de acero; Hormigon 



%%===============/////////===========/////////==============%% 

%% XX Neopreno X///////X X///////X XX %% 

%%%%%%%%%%% %%%%%%%%% V. Acero %%%%%%%%% 



/o /o /o /o /o /o /o /o /o /o /o 



/o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o /o 



% % %%%%%%%%% V. Acero %%%%%%%%% 

% % % % % % % % 

% % % % % % % % 

% % %Pilotes % % % 

% % % % % % 

% ....75m |..12m..|...75m |...12m.|....75m 

% 249 m 

% Estribo Izquierdo Apoyo Apoyo Estribo Derecho 
% (con pilotes) (sin pilotes) 

% 



% Roberto Aguiar 

% CEINCI-ESPE 

% Febrero de 201 0 

% 

% modelo_longitudinal_carrizal_ej1 

% 

% KEI, MEI Matrices de rig y masa de estribo izquierdo (44X44) 

% KEC, MEC Matrices de rig y masa de apoyo central (36X36) 

% KED, MED Matrices de rig y masa de estribo derecho (24X24) 

% KA, MA Matrices de rig y masa de viga de acero (8,8) 

% KT, MT Matrices de rig y masa de estructura sin aporte suelo 

% MATRIZ DE RIGIDEZ Y DE MASA 

[KEI,MEI]=rigidez_estribo_izquierdo_C; % Contribucion de estribo izquierdo 
[KEC,MEC]=rigidez_apoyo_central_C; % Contribucion de apoyo central 
[KED,MED]=rigidez_estribo_derecho_C; % Contribucion de estribo derecho 
[KA,MA]=rigidez_viga_superior_C; %Contribuci6n de viga de acero superior 

% SUBMATRICES CERO 

for i=1 :44; for j=1 :36; CEC(i,j)=0; end; end %Primera fila de Submatrices 

for i=1 :44; for j=1 :24; CED(i,j)=0; end; end 

for i=1 :44; for j=1 :8; CEA(i,j)=0; end; end 

CCEC=zeros(36,36);%Segunda fila de submatrices 

for i=1 :36; for j=1 :24; CCED(i,j)=0; end; end 

for i=1 :36; for j=1 :8; CCA(i,j)=0; end; end 

for i=1 :36; for j=1 :24; CCCED(i,j)=0; end; end%Tercera fila 

for i=1 :36; for j=1 :8; CCCA(i,j)=0; end; end 

for i=1 :24; for j=1 :8; CA(i,j)=0; end; end% Cuarta Fila 
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CAA=zeros(8,8);%Quinta fila 

ngl=164;% NUMERO TOTAL DE GRADOS DE LIBERTAD (COORD. PRINCIPALES) 
% MATRIZ DE RIGIDEZ 



KT=[KEI 


CEC 


CEC 


CED 


CEA 


CEA 


CEA; 


CEC 


KEC 


CCEC 


CCED 


CCA 


CCA 


CCA; 


CEC 


CCEC 


KEC 


CCCED CCCA CCCA CCCA; 


CED' 


CCED' 


CCCED' KED 


CA 


CA 


CA; 


CEA' 


CCA' 


CCCA' 


CA' 


KA 


CAA 


CAA; 


CEA' 


CCA' 


CCCA' 


CA' 


CAA 


KA 


CAA; 


CEA' 


CCA' 


CCA' 


CA' CAA 


CAA 


KA]; 


% 


MATRIZ DE MASAS 








MT=[MEI 


CEC 


CEC 


CED 


CEA 


CEA 


CEA; 


CEC 


MEC 


CCEC 


CCED 


CCA 


CCA 


CCA; 


CEC 


CCEC 


MEC 


CCCED CCCA CCCA CCCA; 


CED' 


CCED' 


CCCED' MED 


CA 


CA 


CA; 


CEA' 


CCA' 


CCCA' 


CA' 


MA 


CAA 


CAA; 


CEA' 


CCA' 


CCCA' 


CA' 


CAA 


MA 


CAA; 


CEA' 


CCA' 


CCCA' 


CA' 


CAA 


CAA 


MA]; 



% CONTRIBUCION DE NEOPRENO, DE DUREZA 70, SHORE 70 

kb=324*4; %rigidez de neopreno de 50 por 50 por 6 cm: Con 3 capas de acero de 2 mm. Ga =7.0 
Kg/cm2 

PR=[37 kb; %Numero de grado de libertad y rigidez horizontal 
132 kb; 141 kb; 142 kb; 149 kb; 150 kb; 157 kb; 158 kb]; 

[Kaisla]=resortes_carrizal(ngl,PR);% el programa sirve para FPS y Neopreno. 

% MATRIZ DE RIGIDEZ CON CONTRIBUCION DE NEOPRENO 

KT=KT+Kaisla; 

% 

% Vectores Jh para sismo horizontal 
for i=1 :ngl; Jh(i)=1 ; end; Jh=Jh'; 

% Caracteristicas del suelo solo de un pilote (Considera tres resortes por nudo) 
%[PROP]=[GDL,FH,FV,Mom,dh,dv,dteta,dhu,dvu,dtetau] Modelo Elastoplastic 

% ESTRIBO IZQUIERDO 

[PRO]=[10 22.0 22.0 0.0 0.60 0.60 0.01.70 1.70 0; %Arcilla Blanda 
9 25.0 25.0 0.0 0.60 0.60 0.0 1 .70 1 .70 0; 

8 25.0 25.0 0.0 0.60 0.60 0.0 1 .70 1 .70 0; 

7 160.0 160.0 0 0.075 0.075 0.0 0.090 0.090 0; %Arena meteorizada 
6 260.0 260.0 0 0.070 0.070 0.0 0.090 0.090 0; 

5 129.0 129.0 0 0.80 0.80 0.0 1.00 1.00 0; %Arcilla dura 

4 350.0 350.0 0 0.80 0.80 0.0 1 .00 1 .00 0; 

3 350.0 350.0 0 0.80 0.80 0.0 1 .00 1 .00 0; 

2 350.0 350.0 0 0.80 0.80 0.0 1 .00 1 .00 0; 

1 350.0 350.0 0 0.80 0.80 0.01.00 1.00 0]; 

[PROPESTI]=genera_gdl_resortes(PRO,2,10); 

% PRIMER APOYO 

[PROP1]=[47 25.0 25.0 0.0 0.60 0.60 0.0 1.70 1.70 0; %Arcilla Blanda 
46 129.0 129.0 0 0.80 0.80 0.0 1.00 1.00 0; %Arcilla dura 

45 350.0 350.0 0 0.80 0.80 0.01.00 1.00 0]; 

% Genera la rigidez de los resortes de los otros pilotes de primer apoyo 
[PROPPI1]=genera_gdl_resortes(PROP1,3,12); 

% SEGUNDO APOYO 

[PROP2]=[83 25.0 25.0 0.0 0.60 0.60 0.0 1.70 1.70 0; %Arcilla Blanda 
82 129.0 129.0 0 0.80 0.80 0.0 1.00 1.00 0; %Arcilla dura 

81 350.0 350.0 0 0.80 0.80 0.01.00 1.00 0]; 

% Genera la rigidez de ios resortes de ios otros piiotes de segundo apoyo 
[PROPPI2]=genera_gdl_resortes(PROP2,3,12); 

% VECTOR CON TODOS LOS RESORTES DEL SUELO 

[PROP]=[PROPESTI; PROPPI1 ; PROPPI2];%COntiene gdi y propiedades de todos ios suelos 

% EMPUJE ACTiVO DEL SUELO EN ESTRiBO IZQUIERDO 

FT=zeros(44,1); 

gama=1 .6; % Peso especifico del suelo 
esuelo=10; %Ancho cooperante del suelo 
for i=1:10;HI(i)=i*1.5;FT(i)=gama*(26.03-HI(i))''2*esuelo/2;end 
Hl(1 1 )=1 6.20;FT(21 )=gama*(26.03-1 6.20)''2*esuelo/2; 

for i=1 :6;HI(1 1 +i)=1 6.20+i*0.99; FT(21 +2*i)=gama*(26.03-HI(1 1 +i))''2*esuelo/2;end 
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HI(18)=22.73;FT(35)=gama*(26.03-22.73)''2*esuelo/2; 

for i=1:3;HI(18+i)=22.73+i*1.1; FT(37+2*i)=gama*(26.03-HI(18+i))''2*esuelo/2;end 
qei=KEI\FT; %Despiazamientos por empuje de suelo 

% EMPUJE ACTIVO DEL SUELO EN ESTRIBO DERECHO 

FD=zeros(24,1); 

HD(1)=1.00; FD(2)=-gama*(1 0.84-1 )''2*esuelo/2; 

for i=1:6;HD(1+i)=1.0+i*0.99; FD(2+i*2)=-gama*(10.84-HD(1+i))''2*esuelo/2; end; 
HD(8)=7.54;FD(1 8)=-gama*(1 0.84-7.54)''2*esuelo/2; 

for i=1:3;HD(8+i)=7.54+i*1.1;FD(18+i*2)=-gama*(10.84-HD(8+i))''2*esuelo/2;end 
qed=KED\FD; %Desplazamientos por empuje de suelo 

% VECTOR DE CONDICIONES INICIALES 

for i=1 :74;qin(i)=0;end; qin=qin';for i=1 :22;qfin(i)=0;end;qfin=qfin'; 
q=[qei;qin;qed;qfin];%Desplazamientos debido a empuje de suelo estatico 

% CONDICIONES INICIALES 

%q=zeros(ngl,1); %Vector de desplazamientos 

qp=q; %Vector de velocidades iniciales 

%qa=zeros(ngl,1); %Vector de despla inic. (Coord. Principales) 

qa=q; 

qap=qa; %Vector de velocidades iniciales (Coord. Principales) 

Y=[qa;qap]; % Vector Y de condiciones iniciales 
zeda=0.05; %Factor de amortiguamiento de estructura 
% Cargo Sismo 
load rec_1s3 
nn=length(ux);tu(1 )=0; 
ugx=0.01*ux; %Acelerograma viene en gals 
for iji=1 :nn-1 
iji 

acel=-ugx(iji);tu(iji+1)=iji*dT; 
for jji=1 :6 

[KTT]=resortes_gdl_esmeraldas(q,KT,PROP); 
[T,phi,omega]=orden_eig(KTT,MT); %Frecuencias de vibracion 
% MATRIZ DE AMORTIGUAMIENTO 
[CT]=amortiguamiento(MT, phi, omega, zeda); 

% Respuesta en el tiempo para sismo horizontal 
[Yn]=pee_de_uno(KTT,CT,MT,Jh,acel,dT,Y); 
qa=Yn(1:ngl);% Nuevos desplazamientos. Coordenadas Principales 
q=qa;% Nuevo vector de desplazamientos 
Q=KT*q;% Vector de Cargas Generalizadas 
Y=Yn; %Se actualiza el vector Y 
end 

for i=1 :ngl 

qt(i,iji+1)=q(i); % Se almacena la respuesta en el tiempo 
Qt(i,iji+1)=Q(i); % Se almacena la respuesta en el tiempo 
end 
end 

% Valores Maximos en absolute 
for i=1 :ngl 

qtmax(i)=max(abs(qt(i, :)));% Desplazamientos maximos 
Qtmax(i)=max(abs(Qt(i, :)));% Cargas maximas 
end 

% DESPLAZAMIENTOS MAXIMOS EN ESTRIBOS 

qmaxEi=[qtmax(1 ) ;qtmax(2) ;qtmax(3) ;qtmax(4) ;qtmax(5) ;qtmax(6) ;qtmax(7) ; 
qtmax(8) ;qtmax(9) ;qtmax(1 0) ;qtmax(21 ) ;qtmax(23) ;qtmax(25) ;qtmax(27) ; 
qtmax(29) ;qtmax(31 ) ;qtmax(33) ;qtmax(35) ;qtmax(39) ; 
qtmax(41 ) ;qtmax(43)] ; 

qmaxEd=[qtmax(1 1 8);qtmax(1 20);qtmax(122);qtmax(1 24);qtmax(1 26); 
qtmax(1 28) ;qtmax(1 30) ;qtmax(1 34) ;qtmax(1 36) ;qtmax(1 38) ; 
qtmax(140)j; 

% FUERZAS MAXIMAS EN ESTRIBOS 

FmaxEi=[Qtmax(1 ) ;Qtmax(2) ;Qtmax(3) ;Qtmax(4) ;Qtmax(5) ;Qtmax(6) ;Qtmax(7) ; 
Qtmax(8) ;Qtmax(9) ;Qtmax(1 0) ;Qtmax(21 ) ;Qtmax(23) ;Qtmax(25) ;Qtmax(27) ; 
Qtmax(29) ;Qtmax(31 ) ;Qtmax(33) ;Qtmax(35) ;Qtmax(39) ; 

Qtmax(41 ) ;Qtmax(43)] ; 

FmaxEd=[Qtmax(1 1 8) ;Qtmax(1 20) ;Qtmax(1 22) ;Qtmax(1 24) ;Qtmax(1 26) ; 
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Qtmax(1 28) ;Qtmax(1 30) ;Qtmax(1 34) ;Qtmax(1 36) ;Qtmax(1 38) ; 

Qtmax(140)]; 

% DESPLAZAMIENTOS MAXIMOS EN APOYOS 

for i=1:12;Z(i)=4.91*i;end 

qmaxCI =[qtmax(45) ;qtmax(46) ;qtmax(47) ;qtmax(48) ;qtmax(49) ;qtmax(50) ;qtmax(51 ) ; 

qtmax(52) ;qtmax(53) ;qtmax(54) ;qtmax(55) ;qtmax(56)] ; 
qmaxC2=[qtmax(81 ) ;qtmax(82) ;qtmax(83) ;qtmax(84) ;qtmax(85) ;qtmax(86) ;qtmax(87) ; 

qtmax(88) ;qtmax(89) ;qtmax(90) ;qtmax(91 ) ;qtmax(92)] ; 

% FUERZAS MAXIMOS EN APOYOS 

FmaxCI =[Qtmax(45) ;Qtmax(46) ;Qtmax(47) ;Qtmax(48) ;Qtmax(49) ;Qtmax(50) ;Qtmax(51 ) ; 

Qtmax(52) ;Qtmax(53) ;Qtmax(54) ;Qtmax(55) ;Qtmax(56)] ; 

FmaxC2=[Qtmax(81 ) ;Qtmax(82) ;Qtmax(83) ;Qtmax(84) ;Qtmax(85) ;Qtmax(86) ;Qtmax(87) ; 
Qtmax(88) ;Qtmax(89) ;Qtmax(90) ;Qtmax(91 ) ;Qtmax(92)] ; 

% 

save TOT 
% Figuras 

figure(1) %Desplazamientos Estribo izquierdo 
subplot(221);plot(tu',qt(10,:)');title('GDL 10'); 
subplot(222) ;plot(tu',qt(37, :)') ;title('G DL 37') ; 
subplot(223);plot(tu',qt(141,:)');title('GDL141'); 
subplot(224) ;plot(tu',qt(1 47, :)') ;title('G DL 1 47') ; 
figure(2) %Desplazamientos Apoyo Central 1 
subplot(221);plot(tu',qt(52,:)');title('GDL 52'); 
subplot(222) ;plot(tu',qt(53, :)') ;title('G DL 53') ; 
subplot(223) ;plot(tu',qt(55, :)') ;title('G DL 55') ; 
subplot(224) ;plot(tu',qt(56, :)') ;title('G DL 56') ; 
figure(3) %Desplazamientos Estribo Derecho 
subplot(221);plot(tu',qt(118,:)');title('GDL118'); 
subplot(222) ;plot(tu',qt(1 32, :)') ;title('G DL 1 32') ; 
subplot(223) ;plot(tu',qt(1 58, :)') ;title('G DL 1 58') ; 
subplot(224) ;plot(tu',qt(1 64, :)') ;title('G DL 1 64') ; 
figure(4) %Fuerzas de Corte Estribo izquierdo 
subplot(221);plot(tu',Qt(10,:)');title('GDL10'); 
subplot(222) ;plot(tu',Qt(37, :)') ;title('G DL 37') ; 
subplot(223);plot(tu',Qt(141,:)');title('GDL141'); 
subplot(224) ;plot(tu',Qt(1 47, :)') ;title('G DL 1 47') ; 
figure(5) %Fuerzas de Corte Apoyo Central 
subplot(221);plot(tu',Qt(52,:)');title('GDL52'); 
subplot(222) ;plot(tu',Qt(53, :)') ;title('G DL 53') ; 
subplot(223) ;plot(tu',Qt(55, :)') ;title('G DL 55') ; 
subplot(224) ;plot(tu',Qt(56, :)') ;title('G DL 56') ; 
figure(6) %Fuerzas de Corte Estribo dERECHO 
subplot(221);plot(tu',Qt(118,:)');title('GDL118'); 
subplot(222);plot(tu',Qt(132,:)');title('GDL132'); 
subplot(223);plot(tu',Qt(158,:)');title('GDL158'); 
subplot(224);plot(tu',Qt(164,:)');title('GDL164'); 
figure(7) % Maximos en Estribos 

subplot(221);plot(qmaxEi',HI);title('Desplazamientos Estribo I.'); 
subplot(222);plot(FmaxEi',HI);title('Fuerzas Estribo I.'); 
subplot(223);plot(qmaxEd',HD);title('Desplazamientos Estribo D.'); 
subpiot(224);plot(FmaxEd',HD);titie('Fuerzas Estribo D.'); 
figure(8) % Maximos en Apoyos Centrales 
subplot(221);plot(qmaxC1',Z);title('Desplazamientos Apoyol.'); 
subplot(222);plot(FmaxC1',Z);title('Fuerzas Apoyo I.'); 
subplot(223);plot(qmaxC2',Z);title('Desplazamientos Apoyo D.'); 
subplot(224);plot(FmaxC2',Z);title('Fuerzas Apoyo D.'); 



% —end— 
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UTILIZACION DE HONGOS MICORRi'ClCO-ARBUSCULARES EN 
PLANTULAS DE TOMATE DE ARBOL (SOLANUM BETACEUM) 
CULTIVADAS IN VITRO: EFECTOS DURANTE LA EASE DE 

ACLIMATACION 

Espm, E., Medina, M.E., Jadan, M. & Proaho, K. 

Carrera de Ingenieria en Biotecnologia 
Escuela Politecnica del Ejercito 



RESUMEN 

El tomate de arbol es una especle de Interes agricola en el pais, que ha demandado en 
los ultimos ahos el desarrollo de nuevas tecnicas de propagaclon. El cultivo in vitro de tejidos 
vegetales constituye una alternativa para la propagaclon del tomate de arbol, ya que ofrece una 
producclon permanente de plantas de calldad y llbres de patogenos. Sin embargo, durante la 
fase de acllmataclon se presenta un alto porcentaje de perdida de las plantas, por su 
Incapacldad de adaptarse a las condiclones amblentales. En este sentido los bongos 
micorricico-arbusculares (HMA), contrlbuyen a mejorar el crecimlento y tolerancia de estas 
plantas, disminuyendo las perdidas durante la fase de acllmataclon. 

En el presente trabajo se apllco un alslado native de HMA en plantulas de tomate de 
arbol cultivadas in vitro, en dos tipos de sustratos: Suelo+Turba y Suelo+Germlnaza. Ademas, 
se estandarizo un metodo de germlnaclon in vitro de semlllas de tomate de arbol, en el que se 
evaluaron: diferentes medlos de cultivo (M&S, M&S con 2mg.L’^ de BAP y M&S con 1,5 mg.L’^ 
de GA 3 ), tipos de semllla (certificadas y de fruto) y tratamlentos pre-germlnativos (escarlficaclon 
mecanica y escarlficaclon mecanica con Imbibiclon). El metodo mas efectivo en la 
estandarizaclon de la germlnaclon in vitro de semlllas de tomate de arbol, fue el tratamlento 
que consistio en semlllas de fruto, escarlficaclon mecanica e Imbibiclon y medio M&S con 
2mg.L^ de BAP. 

Despues de tres meses del proceso de micorrizaclon durante la fase de acllmataclon, 
las plantas fueron evaluadas en: porcentaje de supervivencia, altura, blomasa aerea, peso 
radical, fndice plastocronico (desarrollo), numero de hojas nuevas, area foliar, contenido de 
agua y de fosforo. Los tratamlentos micorrizados fueron significativamente superlores a los no 
micorrizados. Adiclonalmente, se determine la capacldad Infectiva del alslado native de HMA, 
cuyos porcentajes de colonizaclon micorrfcica fueron altos pese al alto contenido nutritive y de 
fosforo de los sustratos. La mejor colonizaclon micorrfcica se obtuvo con el sustrato 
Suelo+Turba contrlbuyendo a un mayor crecimlento de las plantas. 

Esta Investigaclon demuestra que los HMA contrlbuyen de manera positiva a la 
establllzaclon, desarrollo y crecimlento de las plantas de tomate de arbol durante la fase de 
acllmataclon. 



ABSTRACT 

The tree tomato (Tamarlllo) crop Is Important In Ecuador. In recent years. It has required 
the development of new techniques for propagation of this species. The in vitro culture of plant 
tissue Is an alternative to produce tree tomato, as It provides a permanent production of 
pathogen free, quality plants. However, during the acclimatization phase a high number of 
plants die, because of their Inability to adapt to environmental conditions. It has been reported, 
that arbuscular mycorrhizal fungi (AMF), Improve plant growth and tolerance, during the 
acclimatization phase. 

In the present study, we applied a native Isolated of arbuscular mycorrhizal fungi (AMF) 
to tree tomato seedlings in vitro cultured. The Investigation assessed the efficiency of two 
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substrates, Soil+Peat and Soil+Germinaza. In addition, the study inciuded the standardization of 
an in vitro germination method for tree tomato seeds, we evaiuated the treatments: cuiture 
medium (M&S, M&S with 2mg.L'^ BAP and M&S with 1,5 mg.L’^ of GA 3 ), seed type (fruit seeds 
and certified seeds) and pre-germinative treatment (mechanicai scarification and imbibition pius 
mechanicai scarification). The most effective method for the in vitro germination of tree tomato 
seeds, consisted in fruit seeds, mechanicai scarification pius imbibition and M&S with 2mg.L'^ 
BAP medium. 

The triai in the acciimatization phase was maintained for three months and at the end of 
it, the piants were evaiuated by survivai rate, iength, biomass, root weight, piastochron index, 
number of new ieaves, ieaf area, water and phosphorus content. According to the resuits, 
treatments with mycorrhizae showed higher growth and deveiopment, than non mycorrhizae 
treatments. Additionaiiy, we determined the infectivity of the AMF native isoiated, it showed high 
mycorrhizai coionization rates despite the eievated nutrient and phosphorus content of 
substrates. Regarding the substrate, Soii+Peat promoted the mycorhizai coionization and piant 
growth more than Soii+Germinaza. 

The investigation determined that AMF contributes positiveiy to the stabiiization, 
deveiopment and growth of tree tomato piants during the acciimatization phase. 

1. INTRODUCCION 

La importancia economica dei cuitivo de tomate de arboi y ias grandes expectativas 
de exportacion actuaies generan ia necesidad de potenciar tecnicas ecoiogicas de 
mejoramiento y produccion masiva dei fruto, que permitan cumpiir con ios estandares de 
caiidad y sanidad exigidos por ei mercado internacionai (SICA, 2007). 

El cuitivo in vitro de tejidos vegetales se perfila como una herramienta biotecnologica 
eficiente, que contribuira a la produccion a gran escala de tomate de arboi en el pafs. Sin 
embargo, esta tecnica necesita superar las limitaciones que se presentan en la fase de 
aclimatacion de las plantulas mediante la aplicacion de microorganismos beneficiosos que 
forman parte de la microbiota habitual del suelo (Roveda et ai., 2007), y que normalmente son 
eliminados en estos sistemas de propagacion por el uso de sustratos esteriles (Vidal et at., 
1992; Pedraza et ai., 2001). 

El uso de hongos micorrfcico-arbusculares (HMA) durante la fase de aclimatacion de 
plantas micropropagadas, ha despertado un gran interes en Ios ultimos ahos, debido a Ios 
resultados beneficiosos que se han reportado en varias investigaciones en frutales incluido el 
tomate de arboi (Osorio et ai., 2006). Las ventajas inherentes a la simbiosis les permiten a las 
plantas micorrizadas sobrevivir a la etapa de endurecimiento y adaptarse mejor en campo 
(Vega, 1999; Alarcon etat., 2000). 

Los estudios realizados sehalan que las plantas cultivadas in vitro micorrizadas 
enfrentan mejor la fase de aclimatacion (Vega, 1999). La simbiosis mejora en las plantas la 
captacion y asimilacion de nutrientes de baja movilidad como el fosforo. Adicionalmente reduce 
Ios efectos del estres hfdrico a causa de la rafz poco desarrollada y mejora el funcionamiento 
fotosintetico y estomatico. Los FIMA establecen en la planta un sistema de defense frente a 
patogenos, incrementan el crecimiento foliar y radicular (Thacker & Jasrai , 2002) e inducen la 
produccion de hormones en la planta (Alarcon & Ferrera-Cerrato 1999; Pedraza etai., 2001). 

Este trabajo propone establecer una tecnica que implemente el uso de Ios hongos 
micorrfcico-arbusculares para incrementar el desarrollo, adaptabilidad y supervivencia de las 
plantas de tomate de arboi, durante la fase de aclimatacion. Esta tecnologfa contribuira a 
minimizar las perdidas posteriores ai cuitivo in vitro de tomate de arboi y a ofrecer plantas de 
caiidad con el valor agregado de una mayor capacidad adaptativa en campo. 
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2. MATERIALES Y METODOS 

2.1 Material vegetal y microbiologico 

Las semillas se recolectaron en la provincia de Imbabura, Canton Pimampiro, a partir 
de arboles de tomate de arbol de la varledad amarllla. Las esporas nativas de HMA, se alslaron 
de suelos utlllzados en cultivo de tomate de arbol, en el mismo Canton, medlante la tecnica de 
procesamlento de suelo descrita por Gerdeman y Nicholson en 1963, modificada por Herrera et 
at., 2004. 

2.2 Cultivo in vitro de semillas de tomate de arbol 

Para la germlnaclon in vitro de semillas de tomate de arbol, se evaluaron los factores: 
medio de cultivo (M&S, M&S con 2mg.L’^ de BAP y M&S con 1,5 mg.L’^ de GA3), tlpo de 
semllla (certificadas y de fruto) y tratamlento pre-germlnativo (escarlficaclon mecanica y 
escarlficaclon mecanica con Imbibiclon). Las semillas se desinfectaron con cloro al 1,5% por 5 
mlnutos y alcohol al 70% durante 2 mlnutos. La escarlficaclon se realize con un llgero corte con 
bisturi en la testa de la semllla. Los frascos con las semillas se mantuvieron a 25±3-C de 
temperature, 65% de humedad relative, 14 horas luz y ® 2000 luxes, durante 1 mes. 
FInalmente se sembraron 400 semillas con el metodo estandarizado, que consiste en: en 
semillas de fruto, escarlficaclon mecanica e Imbibiclon y medio M&S. 

2.3 Aplicacion de HMA durante la fase de aclimatacion de plantulas de 
tomate de arbol (Solanum betaceum) 

Cuando las plantas presentaron dos hojas verdaderas y una altura de 7 cm, se 
trasplantaron a los sustratos: Suelo+Turba y Suelo+Germlnaza, prevlamente esterlllzados en 
autoclave a 121 -C, 15 lb de preslon, durante 2 horas. Se ahadleron directamente 20 esporas a 
las rafees de las plantulas con ayuda de una micropipeta. Despues de 3 meses (al finallzar la 
fase aclimatacion), las plantas fueron evaluadas en: porcentaje de supervivencia, altura, 
blomasa aerea, peso radical, fndlce plastocronico (0 desarrollo), numero de hojas nuevas, area 
foliar, contenido de agua y de fosforo. Ademas se determine el porcentaje de colonizaclon 
radical medlante la tecnica de tinclon de rafees de Phillips & Hayman, 1 970. 

3. RESULTADOS Y DISCUSION 

La estandarizaclon del metodo para la germlnaclon de semillas, mostro diferenclas 
significativas entre los tratamlentos. El medio de cultivo mas efectivo fue el M&S con 2mg.L'^ de 
BAP, posiblemente debldo a que este bloregulador promueve el crecimlento celular en los 
cotlledones y ellmlna la latencia de las semillas (Fig. 1). Sin embargo, estas semillas 
desarrollaron callos y un desarrollo morfologico anormal, probablemente porque la presencia 
de BAP en el medio Indujo el crecimlento celular formando callos e Inhiblo el crecimlento apical 
del tallo principal, generando el crecimlento de los brotes laterales (Vega y Prehn, 2005). 

En cuanto al tlpo de semllla, las semillas certificadas mostraron muy bajos niveles de 
germlnaclon. Posiblemente porque estas entran en latencia como respuesta adaptativa a los 
tratamlentos prevlos a los que se someten (Koornneef et at., 2002). La escarlficaclon mecanica 
comblnada con la Imbibiclon fue efectiva ya que el corte rompe la barrera mecanica que 
representa la testa y permite la hidrataclon de la semllla durante la Imbibiclon, esto Induce la 
sallda de la radfcula y da paso a la germlnaclon (Manz et at., 2005). 

Los hongos micorrfeleo arbusculares, contrlbuyeron a mejorar el estado fislologico y 
crecimlento de las plantas de tomate de arbol durante la fase de aclimatacion, como se observe 
en la Table 1, las variables evaluadas son significativamente mayores en las plantas 
micorrizadas, que en las no micorrizadas. 
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Numero de semlllas germinadas por tratamiento 
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Figura 1. Numero de semillas germinadas en cada tratamiento del ensayo de germinacion in vitro 6e 
semillas de tomato de arbol (Soianum betaceum), a las 3 semanas de siembra. 



El establecimiento de la simblosis, Incremento la altura, blomasa, numero de hojas 
nuevas y desarrollo de la planta, probablemente debldo a que permite una mayor absorclon y 
asimllaclon de nutrlentes, estimulando tamblen la producclon de fitohormonas. Adiclonalmente 
mejora la actividad fotosintetica, que es deficlente en las plantas que provlenen de cultivo in 
vitro, por el habito heterotrofico que estas adquieren en el medio de cultivo (Read, 1998). 

La contrlbuclon de las micorrizas para mejorar el desarrollo radical es Importante, ya 
que las plantas provenlentes de cultivo in vitro, tienen un sistema radical debll. En el ensayo se 
observe un mayor peso radicular en las plantas micorrizadas, esto se debe a que los HMA 
Incrementan la ramificaclon de la raiz, reducen el tamano de las raices adventiclas y aumentan 
el numero de estas en la planta. Estas modificaclones generan una raiz mas densa y eficlente 
por la presencia del micello externo que penetra los poros del suelo que son Inaccesibles para 
la raiz. (Vega, 1991, Pedraza et al., 2001). 

Se observe que las plantas micorrizadas presentan algunas ventajas en cuanto a la 
absorclon de agua, debldo a que los bongos mejoran la conexion vascular entre tallo y raiz, la 
agregaclon del suelo, la actividad estomatica y la acclon del acldo abscicico, para evltar la 
deshidrataclon de las plantas (Auge, 2004). El contenido de fosforo no fue significativamente 
mayor en las plantas micorrizadas, posiblemente porque los sustratos tenian un alto contenido 
de fosforo. Los HMA, medlante su micello Incrementan la absorclon de fosforo cuando este no 
esta disponible para las plantas. En este case el elemento estaba disponible para las plantas 
micorrizadas y no micorrizadas, por lo que no se observe esta contrlbuclon (Grace etal., 2009). 

Se evidencia ademas en la Tabla 1, que el sustrato Suelo+Turba, favorece tanto el 
crecimlento de la planta come la colonizaclon micorrfcica, potenclando la acclon del bongo en 
variables como numero de bojas nuevas, peso fresco de la blomasa aerea, peso radical y area 
foliar. 

FInalmente, el Inoculo alslado mostro una gran capacldad Infectiva, lo que demuestra el 
exito del use de microorganismos natives del suelo, que desarrollan ventajas adaptativas y de 
compatibllldad con las plantas. 
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Tabla 1 . Variables de crecimiento y desarrollo de las plantas de tomate de arbol, durante la fase de 
aclimatacion, de acuerdo a los tratamientos (con y sin micorrizas, en los diferentes sustratos). 



Variable en planta 


Suelo+Turba 


Suelo+Germinaza 


Micorrizas 


Control 


Micorrizas 


Control 


Altura(cm) 


35.67 a 


24.31 c 


30.39 b 


19.06d 


indice plastocronico 


8.66 a 


6.24 b 


7.76 a 


5.42 b 


Area foliar (cm^) 


214.44 a 


86.67 c 


133.11 b 


63.78 c 


Biomasa aerea (peso fresco g) 


28.20 a 


16.76 c 


19.17b 


5.63 d 


Biomasa aerea (peso seco g) 


2.65 a 


1.62 b 


1.83 b 


1.05 c 


Numero de hojas nuevas 


4.67 a 


4 a 


4.61 a 


2.61 b 


Peso radical (g) 


25.54 a 


9.8 c 


17.09 b 


8.98 c 


Contenido de agua (%) 


90.62 a 


89.89 b 


90.44 a 


89.24 b 


Contenido de fosforo (%) 


0.44 a 


0.38 a 


0.43 a 


0.37 a 


Porcentaie de colonizacion (%) 


52.78 a 


Oc 


37.78 b 


Oc 



4. CONCLUSIONES 



La tecnica estandarizada para la germinacion in vitro de semillas de tomate de arbol 
(Soianum betaceum), que consistio en: medio basico de Murashige y Skoog, semillas 
provenientes de fruto y escarificacion mecanica e imbibicion como tratamientos pre- 
germinativos; promovio la germinacion in vitro de las semillas y permitio obtener plantas viables 
y resistentes. 

Por otro lado, las especies nativas de bongos micorricico-arbusculares aisladas de 
suelos con cultivos de tomate de arbol, de la provincia de Imbabura, mostraron una gran 
capacidad infectiva y tolerancia a altas concentraciones de fosforo, nitrogeno y materia 
organica. 

Los bongos micorricico arbusculares contribuyeron a mejorar la adaptabilidad, 
crecimiento y desarrollo de las plantas de tomate de arbol, durante la fase de aclimatacion. El 
sustrato Suelo+Turba favorecio la colonizacion micorricica y el crecimiento de las plantas de 
tomate de arbol durante la fase de aclimatacion. 

El uso de bongos micorrfcico-arbusculares, para la aclimatacion de plantas cultivadas 
in vitro, constituye una potencial berramienta biotecnologica natural, de bajo costo, que 
permitira minimizar las perdidas en la produccion de plantas a gran escala, ya que mejoran el 
crecimiento y desarrollo de la planta permitiendole enfrentar de mejor manera la fase de 
aclimatacion. 
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RESUMEN 

En la actualidad, la exigencia de los mercados internacionales para la comercializacion 
de productos agricolas sin quimicos, ha promovido el interes en la fabricacion de tecnologias 
que permitan la generacion de herramientas biotecnologicas sustentables. En este contexto, 
los bongos micorricico arbusculares (HMA) han sido probados a manera de biofertilizante con 
gran exito. 

Estos microorganismos que forman parte de la microbiota comun del suelo, influyen en 
el incremento del desarrollo de los cultivos. Ademas, facilitan la recuperacion de suelos 
erosionados y la estabilizacion de espacios contaminados. La presente investigacion tuvo como 
objeto, el estudio de la actividad de los exudados radicales: formononetin (7-hidroxi, 4'metoxi 
isoflavona), tiamina y acido 3-ind6l butirico en la germinacion de esporas de HMA para 
promover la micorrizacion de plantas de avena, maiz y soya y de esta manera mejorar las 
caracteristicas de los inoculos comerciales ya existentes. 

Los resultados obtenidos en este trabajo mostraron que inocular las plantas con HMA 
pre-germinados, acelera el proceso de simbiosis. La micorrizacion incremento 
considerablemente el desarrollo vegetal en relacion con los controles (plantas sin HMA). 
Adicionalmente, se pudo observar que los niveles de produccion de esporas varian 
dependiendo de la especie vegetal, siendo avena y soya los hospederos que mayor poblacion 
de HMA presentaron al final de esta investigacion, lo que permite sugerirlas como plantas utiles 
para la generacion de inoculo micorricico a gran escala. 

ABSTRACT 

Arbuscular mycorrhizal fungi (AMF) are microorganisms constituting the common soil 
microflora, whose role includes the increase of crop yield and recovery and stabilization of 
contaminated soils or eroded. Currently, the requirement of international markets for the 
marketing of agricultural products which form of culture does not involve the application of 
chemicals, has encouraged interest in making new technologies that allow the generation of 
biotechnological tools and friendly sustainable environment and in this context the AMF play a 
key role. 

For the use of biofertilizer HMA way, it is necessary to study all the parameters needed 
to promote their best performance on the ground. In this sense, this research aimed to study the 
activity of root exudates: formononetin (7-hidroxi, 4’-metoxi isoflavona), thiamine and 3-indole- 
butyric acid on the germination of AMF spores and induction of mycorrhizal plants of corn, 
soybeans and oats, to thus increase the functionality of commercial inocula and identify suitable 
plant species for the multiplication of AMF spores in the beds of fertilizer production. 

The results from this study demonstrate that the formononetin and thiamine produced a 
positive effect on the germination of AMF spores in vitro. Also, it was proved that the inoculated 
plants (corn, soybeans and oats) AMF spores pre-germinated, the symbiosis process 
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accelerates, the greater the effect mainly on the biomass plants treated with AMF spores 
without pre-germination. In all cases controls (plants without AMF) had the lowest range of 
growth. It was also noted that production levels of AMF spores, depends on the host and in this 
sense soybeans and oats are plants where sporulation rates are extremely high, which may 
suggest as useful plants for the generation of mycorrhizal inoculum on big scale. 

1. INTRODUCCION 

Actualmente, se estan implementando nuevas formas de cultivo que reemplazan el uso 
de pesticidas con herramientas biotecnologicas que permiten un incremento en la 
productividad, manteniendo el equilibrio biologico del suelo. Es asf, que el uso de hongos 
formadores de micorriza arbuscular (HMA) ha cobrado gran importancia en la agricultura 
mondial en los ultimos ahos (Roman, 2003). 

Los HMA forman parte de la microbiota comun del suelo y juegan un papel fundamental 
en el desarrollo y adaptacion (fitness) de las plantas. Cuando la simbiosis entre el hongo 
micorrfcico y la planta se establece, el microorganismo se vuelve parte integral de su 
hospedero, asegurando el crecimiento de este bajo distintas condiciones ambientales. 

Es asf que, un amplio range de plantas con interes agronomico dependen de las 
micorrizas para su desarrollo, sobre todo cuando se encuentran sembradas en suelos con 
problemas de productividad, ya que los HMA incrementan la resistencia natural de las plantas 
frente a desequilibrios bioticos (patogenos) y abioticos (estres hfdrico, salino o contaminantes) 
(Atkinson, 2009). For ello, pafses como Espaha, Cuba, Brasil, Colombia y Chile estan 
trabajando en la generacion de biofertilizantes elaborados en base a HMA. 

Estos biofertilizantes consisten en sustratos (suelo, turba, vermiculita, arena, entre 
otros) donde se encuentran suspendidas esporas de HMA en concentraciones identificadas 
(100-4000 esporas por gramo) y ocasionalmente compactadas en medios polimerizantes 
(alginato de sodio). Este tipo de tecnicas han resultado innovadoras y facilmente aplicables, sin 
embargo presentan dos graves limitantes: la actividad del inoculo comercial y el tipo de cultivo 
en donde el biofertilizante puede ser aplicado. 

En cuanto a la actividad, el problema consiste en que los HMA al ser almacenados, 
ingresan en un perfodo de latencia caracterfstico de las estructuras fungicas. Este estadfo se 
desencadena debido a que los HMA al ser simbiontes estrictos requieren de la presencia de un 
hospedero vegetal para mantenerse actives. Una vez que las esporas entran en latencia, 
pueden permanecer asf desde varies dfas, meses e incluso hasta su senescencia y muerte 
(Siqueira, 1983). Por lo tanto, cuando los biofertilizantes no se usan inmediatamente, 
disminuyen su efectividad, factor que representa perdida de tiempo y dinero para el agricultor. 

En lo referente al tipo de cultivo, se sabe que los HMA han side empleados con altos 
rendimientos principalmente en frutales y algunas hortalizas mientras que en sembrfos de “cicio 
corto” como cereales, la efectividad de la micorriza disminuye puesto que el tiempo que los 
HMA tardan en realizar la simbiosis espontaneamente es similar al tiempo que la planta 
permanece en el suelo antes de su cosecha. En ese sentido, el empleo de HMA como 
biofertilizantes para cultivos de cosecha temprana resulta poco efectivo, eliminandose por ello 
un amplio sector de consumidores de esta tecnologfa. 

Una manera util de resolver las dificultades antes mencionadas es identificar 
mecanismos que fomenten la actividad de las esporas, para de esta manera impedir que estas 
entren en latencia y asf reducir el tiempo que necesitan para realizar la simbiosis. En ese 
sentido, se conoce que la actividad de las esporas esta aparentemente ligada a la extension de 
su tube germinal, por lo tanto, identificar las sustancias que estimulan la produccion de micelio 
en los HMA resultarfa una herramienta clave para potenciar el efecto de los biofertilizantes en 
campo. 
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Se han realizado investigaciones que demuestran que flavonoides como el 
formononetin (Baptista & Siqueira, 1994), auxinas como el acldo 3-ldol butirico (Kerstin eta!., 
2007), vltamlnas como la tiamlna (Siqueira & Hubbell, 1982), entre otros muchos exudados, 
Inducen la germlnaclon esporas. Identificar las dosis en las que estos compuestos (tiamlna, 
formononetin y acldo 3-lndol butirico) Inducen la germlnaclon de los HMA genera grandes 
ventajas en la producclon de Inoculos comerclales donde estos bongos se encuentran 
presentes, ya que medlante estas sustanclas se lograrian Incrementos en competitividad y 
versatllldad. Adlclonalmente, el estudio de estos exudados se convierte en una herramlenta 
Importante para la comprenslon del cicio simblotico que actualmente presenta bastante 
ambiguedad. 

2. MATERIALES Y METODOS 

En la presente Investigaclon se realizaron dos ensayos, el primero consistio en Induclr 
la germlnaclon de esporas de HMA in vitro medlante el empleo de formononetin, tiamlna y 
acldo 3-lndol butfrico en distintas concentraclones. 

El segundo ensayo probo la funclonalldad de los exudados radicales (formononetin, 
tiamlna y acldo 3-lndol butirico) en contacto con esporas de HMA en presencia de distintos 
hospederos vegetales (soya, avena y maiz). Adlclonalmente se comparo el potenclal de 
micorrizaclon de esporas pre-germlnadas in vitro con respecto a esporas Inoculadas sin 
germlnaclon. 

2.1 Material biologico 

En el primer ensayo, se empleo un coctel de esporas de HMA alsladas sigulendo la 
tecnica de Gerdemann & Nicholson (1963) del paramo de los lllnizas. Se evaluaron tres zonas 
de Intervenclon antropogenica distintas, para de esta manera encontrar esporas con una 
amplla funclonalldad (resistencia a agroquimicos y a estres hfdrico). 

Para el segundo ensayo (macetas trampa) se emplearon tres tipos de tratamlentos, uno 
sin esporas, otro con esporas pre-germlnadas y finalmente uno con esporas sin germlnar, para 
probar la efectividad de los HMA en el desarrollo vegetal y su actividad en presencia de 
reguladores de la elongaclon de tubo germinal. 

Como hospederos vegetales se emplearon plantas de mafz (Zea mays), soya (Glycine 
max) y avena (Avena sativa). 

2.2 inductores de ia germinacion de esporas de HMA 

Como se puede observer en la Table 1 se utlllzaron analogos de exudados radicales 
distrlbuldos en un diseho de bloques completamente al azar (DBCA) en arreglo factorial de 
4x4x3. 



Tabla 1. Tratamlentos establecidos para un factorial con cuatro niveles de formononetin, cuatro niveles de 
tiamlna y tres niveles de acldo 3-indol butirico. 



Tratamientos 


Formononetin x tiamina x acido 3-indol butirico 


Formononetin (pM) 


Tiamina (mg.L-1) 


Acido 3-indol butirico 
(mg.L-1) 


FnTnIBAn 


F1 


F2 


F3 


F4 


T1 


T2 


T3 


T4 


IBA1 


IBA2 


IB A3 




0 


1.5 


2 


2.5 


0 


0.01 


0.1 


1 


0 


1 


1.5 



n es el valor que pueden adquirir los factores en las distintas variables (4 niveles para formononetin y tiamlna y 3 
niveles para IBA), todos ellos Interactuaron en el diseno. 
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3. RESULTADOS 

Al aplicar formononetin y tiamina en los platos de cultivo se incrementa el range de 
germinacion de las esporas de HMA con respecto a la germlnaclon espontanea. La FIgura 1 
muestra que el tratamlento que posee un mayor rendimlento es F4T4IBA1 (93,33%). En este 
tratamlento, se colocaron 2.5pM de formononetin, 1 mg.L’^de tiamlna y 0 mg.L'^de IBA. Al 
finallzar el ensayo, el tratamlento control F1T1IBA1 alcanzo un range de germlnaclon de 20%, 
lo que demuestra que la presencia de flavonoldes (formononetin) y vltamlnas (tiamlna) 
Incrementa el fndice de germlnaclon de esporas de HMA y por tanto Induce su actividad in vitro. 




Figura 1. Porcentaje de germinacion de esporas de HMA inducidas con formononetin, tiamina e IBA, con 

respecto al control tratado con agua. 

Despues de constatar que el tratamlento F4T4IBA1 posee la mayor influencia sobre las 
esporas de HMA en cajas Petri, se corroboro su efectividad en macetas trampa; es decir, en 
presencia de esporas y de hospederos vegetales (avena, mafz y soya). 

Las plantas de avena, mafz y soya tratadas con bongos micorrfcico arbusculares, 
mejoran su rendimlento con respecto a los controles (Fig. 2), observandose que para las 
variables de peso y altura en todas las plantas donde se establece la simbiosis, se incrementa 
la biomasa notoriamente. Ademas, se puede observer que las plantas inoculadas con esporas 
de HMA pre-germinadas alcanzan los mas altos rendimientos, superando a la micorrizacion 
espontanea. 

En cuanto al porcentaje de colonizacion se puede ver que la planta de mayor 
rendimlento es el mafz. Mientras que los fndices mas altos de esporulacion y de produccion de 
micomasa de ectofito estan presentes en avena. Como se puede constatar en la Figura 2, los 
tres hospederos vegetales presentan una alta afinidad por el HMA permitiendo niveles de 
colonizacion intra-radical y esporulacion sumamente altos (principalmente el cultivo de avena). 
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Control Reguladores: plantas tratadas sin HMA regadas con reguladores de la germinacion de esporas. Control 
Agua: plantas tratadas sin HMA regadas con agua. HMA-SGR: hongos micorrfcico arbusculares sin germinar regados 
con reguladores de la germinacion. HMA-SGA: hongos micorrfcico arbusculares sin germinar regados con agua. HMA- 
GR: hongos micorrfcico arbusculares germinados regados con regladores. HMA-GA: hongos micorrfcico arbusculares 
germinados regados con agua. 

Figura 2. Efecto de la utilizacion de esporas de HMA pre-germinadas y reguladores de la germinacion de 
esporas de HMA en las variables: peso, altura, porcentaje de colonizacion, poblacion total de HMA y 

micomasa de ectofito. 



4. DISCUSION Y CONCLUSIONES 



Los hongos micorrfcico-arbusculares juegan un papel sumamente importante en la 
actividad, productividad y estabilidad del suelo, slendo uno de los microorganismos de mayor 
versatllldad entre la microblota. 

Para que se establezca la simblosis entre la planta y el HMA es necesarlo que se 
activen una serle de sehales bloquimicas generalmente coordinadas por exudados extra- 
radlcales y extra-fungicos. 

En el presente se trabajo, al utlllzar Isoflavonoldes (formononetin) y vltamlnas (tiamlna) 
para Induclr la germinacion de esporas de HMA se lograron rendimlentos comprendidos entre 
el 26,67 (F3T3IBA1 y F3T2IBA1) y el 93,33% (F4T4IBA1) mientras que las esporas a las que 
no se apilearon estas sustanclas (germinacion espontanea) muestran un porcentaje de 
germinacion cercano al 20%. Valores simllares se han reportado en varlas Investigaclones 
realizadas por Pereira da Silva & SIqueIra (1997); Akiyama, MatsuzakI & Hayashl.(2005); 
Barreto de Novals (2008), entre otros. 
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La utilizacion de esporas de HMA incrementa el vigor (biomasa) de las plantas de 
avena, soya y maiz (Pereira da Silva & Siqueira, 1997; Azcon Aguilar & Barea, 1982). Cuando 
se emplean esporas de HMA pre-germinadas se evidencia mayor produccion de biomasa 
vegetal, siendo superior a la que presentan plantas inoculadas con HMA sin germinar. 

Probablemente, el incremento del porcentaje de micorrizacion esta relacionado con las 
necesidades nutricionales que presente la planta y a la efectividad del inoculo de HMA utilizado 
(Nair ef a/.,1991; Pereira da Silva & Siqueira, 1997; Siqueira & Hubbell, 1982; Kiriachek et al., 
2009). Zea mays fue la especie con mayor potencial de colonizacion por HMA (94%) 
posiblemente ya que esta planta requiere altas concentraciones de nutrientes para su 
desarrollo y en este sentido, los HMA incrementan su capacidad de absorcion de vitaminas y 
minerales (Smith & Read, 1997; Pacios, 2005). A continuacion se encuentran Avena sativa 
(88,33%) y Glycine max (72,33%), que tambien presentan indices de simbiosis altos, lo que es 
caracteristico de pastes y granos (Pereira da Silva & Siqueira, 1997). 

El hospedero que presento mayor poblacion de esporas por Kg de sustrato al culminar 
este ensayo fue Avena sativa con 9566,67 esporas por Kg de sustrato. Esto se ha visto en 
ensayos anteriores realizados con plantas de avena, trigo y cebada. Las mencionadas 
especies vegetales, son hospederos utiles para la produccion de esporas de HMA debido a la 
alta provision de nutrientes que entregan al micelio externo (Smith & Read, 1997). Glycine max 
posee una capacidad de produccion de 6133,33 esporas por Kg de sustrato, seguida de Zea 
mays alcanza las 3533,33 esporas por Kg de sustrato. Estas dos plantas muestran un 
comportamiento eficiente para la generacion de esporas, sin embargo, no llegan a ser tan 
productivas como Avena sativa, segun los datos obtenidos por el presente ensayo. 
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